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INTRODUCTION 

Among  the  many  water  uses  in  Montana,  irrigation  and  hydropower  are  far 
and  away  the  most  important  in  terms  of  volume  requirements.  At  many  lo- 
cations these  two  uses  compete  for  the  limited  resource,  and  consequently, 
each  is  affected  by  the  activities  of  the  other.  Water  diverted  for  irriga- 
tion is  in  part  consumed  and  thus  not  available  for  power  production.  Even 
the  unconsumed  portion  of  the  irrigation  diversion,  the  return  flow,  is 
delayed  significantly  and  therefore  temporarily  unavailable  for  power  gen- 
eration. Power  production,  although  not  a  consumer  of  water,  does  place  a 
time  demand  on  the  resource  in  that  water  stored  in  reservoirs  for  power 
production  is  temporarily  unavailable  to  irrigators. 

In  this  competitive  setting,  any  change  in  either  user's  activity  almost 
invariably  has  an  impact  on  the  other  user.  Such  changes  are  certain  to  occur 
in  the  future.  For  example,  gradual  and  inevitable  use  changes  will  be 
associated  with  improving  irrigation  efficiencies.  Efficiencies,  which  have 
been  increasing  with  expanded  sprinkler  use  and  reduced  conveyance  loss,  are 
projected  to  increase  even  more  rapidly  in  the  future  (1).  The  increased 
water  supply  derived  from  such  efficiency  improvements  might  conceivably  be 
used  in  a  variety  of  ways  and  thus  have  a  variety  of  impacts.  However  it  is 
used,  it  will  almost  certainly  have  an  impact  on  all  users.  If,  for  instance, 
it  were  used  to  irrigate  additional  lands,  the  same  total  diversion  applied 
now  to  a  larger  area  would  yield  a  reduced  return  flow.  The  reduced  return 
flow  would  negatively  impact  both  downstream  power  producers  and  irrigators. 
In  addition  to  these  efficiency  related  changes  in  water  use,  other  use 
changes  are  likely  to  be  proposed  (as  they  have  been  in  the  past)  for  a 
variety  of  reasons  (economic,  environmental,   technological,  etc.).   These 


proposed  changes,  whatever  their  nature,  would  almost  certainly  impact  users 
to  some  degree. 

There  are  costs  and  benefits,  monetary  and  otherwise,  that  accrue  to 
competing  users  from  any  water  use  or  change  in  water  use.  Competing  water 
users  need  to  understand  the  interdependencies  among  competing  uses  so  that 
they  can  anticipate  the  costs  and  benefits  that  will  result  from  changes  in 
use.  Unfortunately,  there  is  often  a  negative  knee-jerk  reaction  to  any 
proposed  use  or  change  in  use  by  others.  This  sort  of  response  may  not  be 
warranted  since  the  proposed  use  or  change  in  use  may  not  be  disadvantageous 
to  competitive  users.  In  fact  a  user  may  benefit  from  changes  proposed  by  a 
competitive  user.  Power  producers,  for  instance,  may  benefit  from  increases 
in  upstream  irrigation  activity  because  irrigation  diversions  and  the  resul- 
tant delayed  groundwater  return  flows  smooth  streamflow  fluctuations.  Thus 
during  the  winter,  when  power  is  generally  at  a  premium,  streamflows  are 
maintained  at  artificially  high  levels  by  groundwater  return  flows  from  the 
irrigation  activities  of  previous  seasons.  The  groundwater  aquifer's  delaying 
action  serves  the  same  purpose  as  the  power  producer's  surface  reservoir  in 
smoothing  streamflow  variations. 

Water  resource  regulators  and  planners  also  need  to  understand  the 
interdependencies  among  competing  uses  so  that  they  can  properly  manage  the 
resource. 

This  paper  presents  a  methodology  for  quantifying  the  water  use  interde- 
pendencies and  thus  a  mechanism  for  assessing  the  costs  and  benefits  resulting 
from  water  use  changes. 


METHODOLOGY 

The  methodology  was  developed  for  the  Canyoa  Ferry  Reservoir  and  drainage 
basin  with  its  ongoing  irrigation  and  hydropower  operations.  However,  it  is 
quite  general  and  thus  applicable  wherever  agriculture  and  hydropower  compete 
for  water.   Figure  1  is  a  schematic  of  the  Canyon  Ferry  drainage  system. 


Canyon   ^ -v^\  Canyon 

Canyon^^  Ferry 

"        Uuttiow-* 
rvoj  " 


Figure  1.   Schematic  of  the  Canyon  Ferry  Drainage. 


The  irrigation  operation  in  the  basin  is  obviously  dependent  upon  the 
diversion  of  water  to  the  Irrigated  area.  The  hydropower  operation  is  depen- 
dent upon  Canyon  Ferry  Reservoir  inflow.  However,  since  this  inflow  is 
derived  from  the  irrigated  area  return  flow  and  bypass  (streamflow  not  divert- 
ed) ,  the  hydropower  operation  is  in  effect  also  dependent  upon  the  irrigation 
diversion. 


A  methodology  was  developed  to  simulate  the  basin's  irrigation  and 
hydropower  operations  so  that  the  impact  of  changes  in  irrigation  activity 
(changes  in  irrigation  efficiency  or  acreage  irrigated)  might  be  assessed. 
Important  features  of  the  methodology  are  discussed  in  the  following  sections 
with  particular  attention  given  to  the  development  of  the  methodology  and  to 
the  assumptions  employed  in  that  development.  Such  assumptions,  though 
limiting  by  nature,  are  absolutely  necessary  in  a  large  and  diverse  hydrologic 
system  such  as  the  Canyon  Ferry  catchment. 

Simulation  of  the  system  of  Figure  1  is  accomplished  in  the  following 
manner.  The  hydrologic  response  of  the  aggregate  irrigated  area  is  determined 
for  a  specified  sequence  of  irrigation  diversions.  The  irrigated  area  out- 
flows thus  obtained  are  taken  into  the  Canyon  Ferry  Reservoir  where  reservoir 
operation  procedures  and  mass  balance  considerations  determine  their  ultimate 
disposal.  Key  elements  of  the  simulation,  i.e.,  area  aggregation,  the  ir- 
rigated area  response,  the  reservoir  mass  balance  requirement,  and  the  reser- 
voir operation  procedures,  are  discussed  under  separate  headings  below. 

AGGREGATION  OF  IRRIGATED  AREAS 

As  indicated  above,  and  as  suggested  by  Figure  1,  the  irrigated  areas 
within  the  Canyon  Ferry  drainage  are  considered  in  aggregate.  Although  this 
idealization  effectively  ignores  the  irrigation  caused  spatial  and  temporal 
variations  in  water  distribution  that  actually  exist  within  the  basin,  it  does 
not  materially  affect  the  hydrologic  behavior  of  the  basin  as  a  whole  (the 
concern  of  this  study).  This  fact  can  be  seen  from  a  consideration  of  Figure 
2. 


Figure  2.   The  effect  of  aggregating  irrigated  areas. 


An  Irrigation  diversion  at  A  reduces  the  streaoflow  at  all  points  below  A. 
The  subsequent  return  flow  from  that  diversion  increases  the  flow  at  all 
points  below  A.  If  the  same  diversion  and  irrigated  area  were  at  B  (instead 
of  A) ,  the  streamflow  impact  would  be  the  same  but  it  would  be  felt  only  at 
points  below  B.  Thus,  the  streamflow  within  the  basin  (between  A  and  B)  is 
materially  affected  by  the  change  in  the  point  of  diversioQ.  However,  the 
streamflow  from  the  basin  as  a  whole  (at  C)  is  affected  only  to  the  extent 


that  the  impact  of  the  diversion  and  return  flow  is  felt  one  day  later  when 
the  diversion  is  made  at  the  same  point  in  time  but  at  A  rather  than  B.  This 
is  due  to  the  greater  instream  travel  time  from  A  (Figure  2). 

Within  the  Canyon  Ferry  catchment,  instream  travel  times  from  the  various 
important  irrigation  sites  are  quite  similar.  It  is  estimated  that  instream 
travel  times  to  Canyon  Ferry  from  most  of  the  important  irrigation  sites  are 
on  the  order  of  1  to  4  days  (a  distance  of  50  miles  per  day  is  traversed  with 
a  stream  velocity  of  3  ft/s).  Thus,  placement  of  the  aggregate  irrigated  area 
at  a  location  about  2  days  distant  from  Canyon  Ferry  would  result  in  flow 
response  discrepancies  of  from  1  to  2  days  for  some  irrigated  areas. 

The  effect  of  aggregation  and  the  resultant  streamflow  response  discrep- 
ancy must  also  be  viewed  in  terms  of  the  time  frame  selected  for  simulation. 
A  response  discrepancy  of  1  to  2  days  Is  very  significant  if  the  simulation 
interval  is  a  day  or  even  a  week.  It  is  much  less  significant  if  the  interval 
is  a  month.  In  this  study,  a  time  interval  of  one  month  is  used  and  the 
travel  time  from  the  aggregate  irrigated  area  is  assumed  to  be  zero  days, 
i.e.,  water  is  assumed  to  flow  directly  from  the  aggregate  area  into  Canyon 
Ferry  Reservoir. 

The  aggregation  of  Irrigated  areas  and  neglect  of  instream  travel  time 
greatly  simplify  the  analysis  with  minimal  loss  in  rigor.   Furthermore,  they 
are  dictated  by  data  inadequacies.   The  daily  data  (streamflow,  etc.)  required 
to  properly  assess  instream  travel  time  delays  between  irrigated  areas  are  not 
available. 

HYDROLOGIC  RESPONSE  OF  THE  IRRIGATED  AREA 

Simulation  of  the  hydrologic  response  was  accomplished  by  requiring  that 


the  continuity  equation  and  the  return  flow  equation  be  satisfied  for  the 
aggregate  irrigated  area.   The  development  and  utilization  of  these  equations 
are  described  in  the  following  sections. 
Continuity  Equation 

Mass  (water)  is  conserved  in  the  hydrologic  system  of  Figure  1.  The 
continuity  equation,  a  mathematical  statement  of  this  fundamental  physical 
principle,  thus  provides  a  basis  for  characterizing  the  systems  hydraulic 
behavior.  The  continuity  or  balance  equation  written  for  the  reach  of  streaa 
traversing  the  aggregate  irrigated  area  (Figure  3)  is 


Groundwater- 
Natural  (GN) 


Tributary 
Inflow  (TI) 


Channel 
Precipitation  (P) 


Inflow  {I)  -¥ 
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Channel   t  Storage  (S)  ^ 


'Outflow 


Diversion  (D)  Surface  Runoff- 

Irrigation  (SRI) 
Groundwater- 
Irrigation  (GI) 


Figure  3.   Stream  reach  through  aggregate  irrigated  area. 


i  +  P  +  Ti  +  SRN  +  sfti  +  Gn  +  GI  -  6  -  6  -  ^ 

where  I.  P,  etc.  are  instantaneous  volumetric  flow  rates.  Note  that  surface 
runoff  and  groundwater  flow  from  natural  processes  are  identified  separately 
from  those  resulting  from  irrigation  activity.  The  equation  states  that  the 
total  instantaneous  volumetric  inflow  rate  to  the  reach  minus  the  total 
instantaneous  volumetric  outflow  rate  from  the  reach  must  equal  the 
instantaneous  rate  of  change  of  water  storage  in  the  reach. 

Because  instantaneous  rates  are  not  a  concern  in  this  study  (nor  are  such 
measurements  generally  available),  the  equation  is  approximated  as 

I  +  P  +  TI  +  SRN  +  SRI  +  GN  +  GI  -  D  -  0  =  -^ 

where  I.  P.  etc.  are  time-averaged  flow  rates  over  the  interval  At.   The 

equation  may  also  be  written 

fl  +  P  +  TI  +  SRN  +  SRI  +  GN  +  GI  -  D  -  0)At  =  AS  =  S2-S, 
where  S^  and  S^  represent  respectively  the  water  volume  stored  in  the  channel 

reach  at   the  beginning  and  end  of   the   time   interval  At.   When  the 

time-averaged  volumetric  flow  rates  within  the  parentheses  are  multiplied  by 

At,  a  volumetric  relationship  results. 

I  +  P  +  TI  +  SRN  +  SRI  +  GN  +  GI  -  D  -  0  =  S^-S^ 
where  I.  P.  etc.  are  volumetric  accumulations  during  At.   It  is  this 
volumetric  form  of  the  continuity  or  balance  equation  that  is  employed  in  this 
study.  As  previously  indicated,  the  time  interval  At  is  a  month.   Thus  I.  for 
example,  represents  the  volumetric  inflow  to  the  reach  over  a  month's  time. 

The  equation  may  be  used  to  determine  any  one  term  provided  all  of  the 
other  terms  are  known.   Unfortunately  numerous  terms  in  the  equation  are 


unknown.   Fortunately  not  all  terms  are  of  equal  significance  (magnitude)  and 
consequently  their  neglect  does  not  materially  disturb  the  equality. 

Because  channel  storage  within  the  irrigated  reach,  S,  is  not  known  nor 
easily  determined,  the  change  in  storage,  S_-S.,  on  the  right  side  of  the 
equation  is  neglected.  The  significance  of  this  neglect  can  be  qualitatively 
assessed  as  follows.  If  the  travel  time  through  the  somewhat  nebulous  stream 
reach  of  Figure  3  is  on  the  order  of  1.5  days,  then  the  reach  is  refilled 
approximately  20  times  each  month  and  inflow  is  20  times  larger  than  storage. 
The  change  in  channel  storage,  S»-S.,  is  of  course  considerably  less  than  S. 
Thus  it  can  be  argued  that  the  storage  change,  S„-S.,  on  the  right  side  of  the 
equation  is  relatively  quite  small  compared  to  the  inflow  on  the  left  side  and 
can  therefore  be  neglected  without  greatly  disturbing  the  equality. 

Similar  relative  magnitude  arguments  could  be  presented  to  justify  the 
neglect  of  channel  precipitation  P,  natural  surface  runoff  SRN,  and  natural 
groundwater  flow  GN  within  the  limited  reach  of  Figure  3.  Fortunately,  as 
will  be  demonstrated  later,  it  is  not  necessary  to  ignore  these  quantities 
even  though  individually  they  are  unknown.   The  continuity  equation  is  now 

(I  +  TI  +  P  +  SRN  +  GN)  +  (SRI  +  GI)  -  D  -  0  =■  0  (1) 

The  terms  are  rearranged  and  the  parentheses  introduced  to  facilitate  further 
discussion. 

The  terms  within  the  first  set  of  parentheses  represent  natural  inflow  to 
the  reach  while  the  terms  in  the  second  set  represent  irrigation  caused 
inflow.  Irrigation  return  flow  (SRI  +  GI)  is  derived  from  and  dependent  upon 
the  irrigation  diversion  D.  The  mathematical  characterization  of  this 
physical  dependency,  an  essential  element  of  the  hydrologic  simulation 
methodology,  is  discussed  in  the  following  section. 


Irrigation  Return  Flow 

The  application  of  irrigation  diversion  water  is  not  100  percent 
efficient.  A  portion  of  the  applied  water  typically  runs  off  of  the  surface 
while  another  portion  percolates  beneath  the  root  zone  to  recharge  the 
groundwater  system.  These  waters  ultimately  reach  the  stream  as  the 
previously  identified  irrigation  surface  runoff  SRI  and  irrigation  groundwater 
flow  GI  in  the  continuity  equation.  The  magnitude  of  the  groundwater  and 
surface  water  irrigation  return  is  dependent  upon  the  magnitude  of  the 
irrigation  diversion  and  upon  the  efficiency  and  method  of  application  of  the 
diversion.  Proper  evaluation  of  the  continuity  equation  is  dependent  not  only 
upon  the  magnitude  of  these  two  return  flows  but  also  upon  the  timing  of  the 
returns.  It  is  the  latter  of  these  dependencies,  the  timing  of  the  returns, 
that  is  addressed  in  the  following  paragraphs. 

The  timing  of  surface  returns  from  irrigation  applications  is  easily 
dealt  with.  It  is  assumed  that  the  overland  flow  time  from  the  point  of 
irrigation  water  application  to  the  stream  is  negligible.  Thus  the  surface 
return  from  a  given  diversion  appears  at  the  stream  during  the  same  month  in 
which  the  diversion  is  made.  In  view  of  the  close  proximity  of  the  irrigated 
areas  to  the  streams  in  the  Canyon  Ferry  drainage,  this  assumption  is  believed 
to  be  quite  reasonable  and  thus  a  very  minor  source  of  simulation  error. 

The  treatment  of  groundwater  returns  is  more  complex  because  the  returns 
are  much  slower  and  continue  over  a  much  longer  period  of  time.  The 
groundwater  return  from  a  given  irrigation  application  decreases  exponentially 
with  time  and  thus  supplements  streamflow  for  many  months  after  the  month  of 
irrigation  application. 
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Glover  (2)  developed  a  procedure  for  decemlnlng  groundwater  return  that 
Is  based  on  a  parallel  drain  concept.  Hurley  (3)  used  the  procedure  to 
qxiantlfy  return  flow  to  the  Rio  Grande  River  from  the  Irrigated  Mesllla  Valley 
In  New  Mexico.  The  procedure  yielded  a  return  flow  pattern  that  agreed  quite 
well  with  the  observed  return  flow  pattern.  The  physical  basis  for  Glover's 
procedure  Is  presented  In  the  following  paragraphs. 

When  a  differential  form  of  the  one-dlmenslonal  continuity  equation  Is 
combined  with  the  Dupult-Forchheimer  assumption  for  the  control  volume  between 
the  parallel  drains  of  the  Isotropic,  homogeneous  aquifer  of  Figure  4.  the 
resulting  equation  may  be  written 


Water  Table  Immediately 
After  Irrigation  (at  E=o) 

ater  Table  at  time  t 


Impermeable 


Figure  4.   Parallel  drains, 
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K   3_  ..  ,  dh^.      ^        8h' 
S   3x  ^    8x  ^  dt 


S  Ix  ^^-'^^  3^  =  3F 

where  K  is  aquifer  hydraulic  conductivity  and  S  is  aquifer  storage  coefficient 
(specific  yield).  If  h,  the  drainable  depth,  is  small  compared  to  D,  as  it  is 
in  many  field  situations,  the  nonlinear  equation  can  be  linearized  by  neglect- 
ing h 

K  D   32h  ^  3h  (2) 

S  TP  "   3t 
For  the  conditions 

h  =  0  when  x  =»  0  for  t  >  0 

h  =  0  when  x  =  L  for  t  >  0 

h  =  H  for  0  <  x  <  L  when  t  =  0 

the  solution  for  h  at  any  time  t  and  at  any  location  x  is 


niT 
sin  (—  x)  (3) 


KD   T 
where  a  =  —  =  —  .   Aquifer  transmissivity  is  denoted  by  T. 
•^         S 


This  solution  yields  the  water  table  position  h  at  any  time  after  an  irriga- 
tion application  that  instantaneously  raises  the  water  table  a  distance  H 
(Figure  4) . 

Now  suppose  that  instead  of  an  instantaneous  recharge,  a  continuous 
recharge  at  the  rate  i  occurs.  The  recharge  rate  i  is  the  rate  at  which  water 
is  reaching  the  water  table.  The  units  of  i  are  length  divided  by  time 
(inches  per  month  for  example).  If  t  is  a  time  variable  running  from  o  to  t , 
then  the  increment  of  recharge  (water  table  rise  dH)  occurring  at  any  time  t 


\ 

. 

2   2 

n  TT  at 

n 

h  =  AH; 

^     n=l,3 

5.. 

e 

may  be  written  as  (i/s)dT.   Paralleling  the  development  of  Equation  (3),  the 
drainable  depth  dh  at  time  t  from  this  recharge  increment  would  be 


00 

n-ir-  a(t-T) 

\ 
/ 

e 

L^ 

n-1,3,5.. 

n 

dh  =  ^  i      Z_        ^ sin  (^x)dT 

TT   S 
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For  a  constant  recharge  rate  i,  integration  yields 
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for  the  water  table  position  h  at  any  time  t  and  at  any  location  x. 

By  the  Dupuit-Forchheimer  assumption,  the  flow  to  a  unit  length  of  drain 
(from  both  sides)  is 

Q=2KD   (|^)_^^^  (5) 


^th 
The  water  table  gradient,  (^— )   -,   required  to  determine  drain  flow  mav  be 

°  'x  x=0 

obtained  by  differentiating  Equation  (4)  with  respect  to  x  and  evaluating  the 
derivative  at  x  =  o.  Substitution  of  the  gradient  so  obtained  into  equation 
(5)  yields  the  inflow  to  a  unit  length  of  drain  at  any  time  t, 
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The  cumulative  inflow  volume  per  unit  length  of  drain  at  any  time  t  after  the 
continuous  recharge  starts  is  obtained  by  integrating  Q  with  respect  to  t  from 
o  to  t. 


/  Qdt 
Dividing  the  cumulative  volume  of  drain  inflow  at  time  t  by  iLt,  the  cumula- 
tive groundwater  recharge  at  time  t,  yields 
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Inspection  of  this  equation  reveals  that  R  is  a  function  of  a,  t  and  L.  The 
equation  has  been  evaluated  and  the  results  tabulated  in  terms  of  ^(2),  The 
tabulation  provides  a  means  for  determining  drain  inflow  (return  flow)  when 
groundwater  recharge  is  continuous  at  the  constant  rate  i.  As  can  be  seen 
from  Equation  (6),  an  R  value  (from  the  tabulation)  represents  the  cumulative 
return  flow  volume  at  time  t  expressed  as  a  percentage  of  the  cumulative 
recharge  at  time  t. 

Equation  (6)  does  not  directly  yield  return  flows  for  the  more  common 
situation  when  recharge  is  neither  continuous  nor  at  a  constant  rate.  However 
it  does  provide  the  basis  for  determining  such  returns.  The  principle  of 
superposition,  which  states  that  solutions  to  linear  equations  are  additive, 
is  employed  in  dealing  with  this  situation.  Consider  the  continuous 
groundwater  recharge  of  Figure  5. 
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Figure  5.   Continuous  groundwater  recharge  at  a  constant  rate  1. 

The  recharge  period  may  be  divided  into  equal  intervals  (periods  1,2,3, 

in  this  example).  Because  Equation  (2)  Is  linear,  return  flows  from  the 
separate  periods  are  additive.  Stated  another  way,  the  return  flow  from  a 
particular  period  is  independent  of  the  previous  return  flow  history.  It 
follows  that  a  pare  of  the  return  flow  that  appears  during  period  3  is  from 
the  recharge  during  period  1,  part  is  from  the  recharge  during  period  2,  and 
part  is  from  the  recharge  during  period  3.  If  the  portion  of  the  period  1 
recharge  that  appears  during  period  j  is  designated  r(i,j),  the  cumulative 
return  flow  at  the  end  of  the  third  period,  R.,  can  be  written 

R^  «  rd.l)  +  r(1.2)  +  r(l,3)  +  r(2,2)  +  r(2,3)  +  r(3.3)    (7) 
Similarly,  the  cumulative  return  flow  at  the  end  of  period  2  is 

R2  -  r(l,l)  +  r(1.2)  +  r(2,2)  (8) 

and  at  the  end  of  period  1  is 

R^  -  r(l.l)  (9) 

The  values  R.,  R-,  and  R-  can  be  determined  from  Equation  (6).  A  manipulation 
of  these  three  equations  yields  the  value  r(l,3),  the  return  flow  during 
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period  3  from  the  recharge  of  period  1.  As  will  be  shown,  this  and  similarly 
determined  quantities  r(l,2),  r(l,4),  r(l,5),  etc.  are  essential  elements  in 
the  determination  of  return  flows  from  groundwater  recharge  that  is  neither 
continuous  nor  at  a  constant  rate. 

The  value  of  r(l,3)  can  be  obtained  in  the  following  manner.  Subtraction 
of  the  cumulative  return  flow  of  Equation  (8)  from  the  cumulative  return  of 
Equation  (7)  yields 

R3-R2  -  AR3_2  -  r(l,3)  +  r(2,3)  +  r(3,3)  (10) 

The  quantity  AR_  „  is  simply  the  total  return  flow  that  appears  during  the 
third  period.  It  consists  of  return  flows  from  recharges  during  periods  1,2, 
and  3.  Subtraction  of  Equation  (9)  from  Equation  (8)  yields  the  return  flow 
that  appears  during  period  2. 

R2  -  R^  »  AR2_j  -  r(l,2)+r(2,2)  (11) 

It  consists  of  return  flows  from  recharges  during  periods  1  and  2.  As  previ- 
ously pointed  out,  return  flow  response  from  a  particular  recharge  event  is 
Independent  of  the  previous  return  flow  history.  Thus  the  identical  recharge 
events  of  Figure  5  result  in  identical  return  responses,  but  displaced  in 
time.   It  then  follows  that 

r(l,2)  +  r(2,2)  =  r(2,3)  +  r(3,3) 
and  that  Equation  (10)  may  be  written 

AR3_2  -  r(l,3)  +  r(l,2)  +  r(2,2)  (12) 

Subtraction  of  Equation  (11)  from  Equation  (12)  yields 

AR3_2  -  AR2_^  -  r(l,3) 
This  is  the  return  flow  that  appears  during  period  3  from  the  recharge  during 
period  1.   Similar  manipulation  of  these  and  other  R  values  (determined  from 
Equation  (6)  evaluated  at  4t,  5t,  etc.)  yields  additional  r(l,j)  values.   The 
return  flow  volumes  r(l,l),  r(l,2),  r(l,3), ,  thus  obtained  characterize 
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the  return  flow  response  from  the  period  1  recharge  on  a  period  by  period 
basis. 

If  r(l,j)  values  are  determined  for  a  first  period  unit  volume  recharge 
(ILt  -  1),  then  each  r(l,J)  value  can  be  thought  of  as  a  percentage.  That  is, 
r(l,J)  is  the  percentage  of  the  unit  recharge  during  period  1  that  appears  in 
period  j.  When  tabulated,  these  percentages  (factors)  provide  a  means  for 
determining  the  return  flows  when  groundwater  recharge  is  neither  continuous 
nor  at  a  constant  rate.  Consider  the  time  distribution  of  recharge  shown  in 
Figure  6. 


[YJ  2[T]a  ,  5  6  7  rrl  qI  I  \ 


Figure  6.  Noncontinuous  groundwater  recharge  at  a  variable  rate  1. 

If  V  is  the  recharge  volume  during  period  1,  then  the  return  flow  during 
period  j  resulting  from  the  period  1  recharge  is  r(i,j)V  .  The  total  return 
flow  during  period  5,  for  example,  is 

r(l,5)Vj  +  r(2,5)V2  +  r(3.5)V^  +  r(4,5)V^  +  r(5,5)V^ 
As  previously  described,  return  flow  responses  to  identical  recharge  events 
are  Identical  but  displaced  In  time.   Thus 
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r(2,5)  =  r(l,4) 

r(3,5)  =  r(1.3) 

r(4.5)  -  r(l,2) 

r(5,5)  -  rd.l) 
and  the  total  return  during  period  5  is 

r(l,5)Vj  +  r(l,4)V2  +  r(l,3)V3  +  r(l,2)V^  +  rCl.DV^ 
In  the  example  of  Figure  6,  there  was  no  recharge  during  period  4  or  5  and  so 
V,  »  V  ■  0  and  the  total  return  during  period  5  is 

r(l,5)V^  +  r(1.4)V2  +  r(l,3)V2 
From  the  table  of  previously  determined  r(l,j)  values  and  specified  period 
recharge  volumes,  V.,  the  total  return  flow  during  the  5th  period  can  be 
determined.   The  return  flow  during  any  other  period  can  be  determined  in  a 
like  manner. 

Computed  r(l,j)  values,  the  basis  for  the  return  flow  determination, 
become  smaller  as  j  becomes  larger,  but  they  do  not  mathematically  go  to  zero. 
Thus,  return  flows  from  a  particular  recharge  event  become  progressively 
smaller  with  time  but  continue  indefinitely.  From  a  practical  standpoint 
however,  when  r(l,j)  values  and  associated  return  flows  become  very  small,  and 
substantially  all  (say  99.9%)  of  the  event  recharge  has  returned,  it  is 
reasonable  to  assume  that  subsequent  r(l,j)  values  are  zero  and  that  return 
flow  from  that  event  ceases. 

The  return  flow  determination  procedure  just  described  can  be  extended  to 
situations  where  the  recharge  rate  is  continuously  varying  if  the  continuously 
varying  rate  can  be  reasonably  approximated  by  a  stepwise  rate  variation  as 
shown  in  Figure  7  below. 
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Figure  7.   Stepwise  approximation  of  continuously  varying  recharge 
events. 

The  stepwise  approximation  is  of  course  like  the  recharge  variation  of  Figure 
6  and  can  be  dealt  with  In  the  same  way. 

Assuming  that  the  recharge  of  Figure  7  is  the  result  of  an  irrigation 
application  that  is  similarly  distributed  in  time,  it  follows  that  the  previ- 
ously described  computational  procedures  provide  a  link  between  irrigation 
diversions  (applications)  and  groundwater  return  flows.  The  procedures  are 
appropriate  for  an  irrigated  area  bordered  by  parallel  drains  that  are  spaced 
a  distance  L  apart  as  In  Figure  4.  An  equivalent  problem  is  one  in  which  a 
valley  of  width  L  is  drained  by  a  stream  midway  between  the  valley  borders 
(Figure  8) .      - 
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Figure  8.  Parallel  drains  and  stream  valley. 

The  problems  are  equivaleQC  because  in  each  case  the  flow  Is  zero  through  a 
vertical  plane  located  a  distance  L/2  from  the  drain. 

As  will  be  discussed  in  more  detail  later,  the  two  idealized  flow  sit- 
uations of  Figure  8,  along  with  Equation  6  and  the  relationships  derived 
therefrom,  provide  the  basis  for  computing  ground  water  return  flows  in  the 
Canyon  Ferry  drainage. 

Combined  Use  of  the  Continuity  and  Return  Flow  Relationships 

The  continuity  equation,  as  previously  developed  for  the  stream  reach 
traversing  the  aggregate  irrigated  area,  is 

(I  +  TI  +  P  +  SRN  +  GN)  +  (SRI  +  GI)  -  D  -  0  -  0  (1) 
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where  the  terms  are  as  defined  in  Figure  3.  The  equation  may  be  solved  for 
any  one  unknown  provided  all  other  quantities  are  known.  Certain  of  the  terms 
In  the  equation  are  linked  through  the  previously  discussed  surface  and  ground 
water  return  flow  relationships.  Thus  surface  return  from  irrigation,  SRI, 
and  ground  water  return  from  irrigation,  GI,  can  be  detenniaed  if  irrigation 
diversions,  D,  are  known.  It  follows  that  if  total  reach  inflow  (I  +  TI  +  P  + 
SRN  +  GN)  and  diversions  are  known,  then  reach  outflow,  0,  can  be  determined. 

With  reference  to  Figure  3,  it  can  be  seen  that  total  reach  inflow  (I  + 
TI  +  P  +  SRN  +  GN)  consists  of  reach  inflow  I,  tributary  Inflow  TI,  pre- 
cipitation that  falls  directly  on  the  channel  P,  natural  surface  runoff  SRN, 
and  natural  ground  water  recharge  (or  discharge)  GN.  To  directly  determine 
total  inflow  would  require  many  measurements.  Reach  inflow,  for  instance, 
would  require  streamflow  measurements  on  all  streams  that  enter  the  various 
Irrigated  areas.  Natural  ground  water  recharge  (or  discharge)  would  require 
sophisticated  and  difficult  measurements  at  many  locations.  Similar  problems 
would  be  encountered  in  the  direct  determination  of  the  other  components  of 
the  total  inflow.  Thus  direct  determination  of  separate  inflow  components  is 
not  feasible.  However,  indirect  determination  of  the  combined  or  total  inflow 
(the  concern  of  this  study)  is  feasible  via  the  continuity  equation.  As 
previously  stated,  the  continuity  equation  may  be  solved  for  any  term  provided 
all  other  terms  are  known.  Thus  if  reach  outflows  0  and  diversions  D  (and 
therefore  irrigation  returns  SRI  and  GI)  are  known,  total  inflow  can  be 
determined.  In  the  Canyon  Ferry  system,  historical  diversions  and  irrigated 
reach  outflows  are  known  and  consequently  total  inflow  can  be  determined.  The 
inflow  thus  determined  is  the  "natural"  Inflow  to  the  irrigated  area.  Once 
established,  this  "historical"  record  of  inflow  can  be  used  to  assess  the 
Impacts  of  changes  in  irrigation  activity.   This  is  accomplished  by  solving 


the  continuity  equation  again,  but  this  time  for  outflow  0.  In  this  applica- 
tion of  the  continuity  equation,  total  inflows  (computed  total  inflows)  and 
diversions  (real  or  assumed)  are  known.  A  comparison  of  results  for  various 
assumed  irrigation  diversions  and  efficiencies  reveals  the  impacts  of  changes 
in  irrigation  activity. 

Characterization  of  the  irrigated  area's  hydrologic  behavior,  as 
described  in  this  section,  has  necessarily  utilized  many  assumptions.  Where 
appropriate  the  assumptions  have  been  discussed  with  regard  to  their 
reasonableness  and  their  ramifications.  The  reasonableness  of  some  of  the 
assumpClons,  however,  can  only  be  judged  by  comparing  the  assumed  condition  to 
real  physical  condition.  For  example,  the  assumption  of  an  isotropic, 
homogeneous  aquifer  may  be  reasonable  for  some  real  aquifers  and  unreasonable 
for  others.  Until  the  physical  character  of  the  real  aquifer  is  examined,  the 
reasonableness  of  the  assumption  cannot  be  judged.  Further  discussion  of  such 
assumptions  must  therefore  await  the  more  detailed  description  of  the  Canyon 
Ferry  system  presented  in  a  later  section. 

RESERVOIR  MASS  BALANCE 

The  irrigated  area  outflow  of  the  previous  section  flows  into  the  Canyon 
Ferry  Reservoir.   As  with  the  irrigated  area  stream  reach,  mass  (water)  is 
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Figure  9.   Canyon  Ferry  Reservoir. 

conserved  in  the  reservoir  (Figure  9)  and  the  continuity  equation  may  be  written 

RI  +  RP  -  RE  -  RG  -  RO  -  RSj  +  BS^  -  RS.  -  BS.  ' 
where  the  terms  are  as  Identified  in  the  figure.  The  terms  RI,  RP,  etc.  on 
the  left  side  of  the  equation  are  volumetric  accumulations  during  the  month. 
The  terms  on  the  right  side,  RS.,  etc.  are  storage  volumes  at  the  beginning 
and  end  of  the  month.  The  subscript  1  identifies  the  beglnnlng-of-the-month 
storage  while  the  subscript  2  identifies  the  end-of-month  storage.  System 
storage  consists  of  two  components:  reservoir  storage  RS  and  bank  storage  BS. 

The  equation  can  be  solved  for  any  one  term  provided  all  other  terms  are 
known.  Thus  a  rigorous  solution  demands  that  monthly  values  be  known  for  all 
of  the  independent  variables  for  the  period  of  Interest.  Because  such  a  data 
base  does  not  exist  for  the  Canyon  Ferry  system,  some  rigor  must  be 
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sacrificed.  The  following  paragraphs  identify  the  terms  in  the  continuity 
equation  that  were  neglected  in  this  study.  The  significance  of  that  neglect 
is  also  addressed. 

A  complete  lack  of  information  about  bank  storage  made  it  necessary  to 
neglect  bank  storage  change  (BS-  -  BS.).  Bank  storage  changes  with  reservoir 
level.  A  rise  in  reservoir  level  produces  an  hydraulic  gradient  that  causes 
water  to  flow  into  the  adjacent  soil  where  it  is  stored  until  a  falling 
reservoir  level  reverses  the  gradient  and  causes  water  to  flow  from  the  soil 
into  the  reservoir.  The  amount  of  water  stored  in  the  soil  depends  upon  the 
reservoir  level,  the  soil's  storage  coefficient  (specific  yield),  and  the 
extent  of  the  soil  as  determined  by  the  bed  rock  configuration.  Although  the 
volume  stored  may  be  quite  large,  the  monthly  change  in  volume  stored  (BS„  - 
BS.)  is  much  less.  The  monthly  volume  change  is  dependent  upon  the  monthly 
reservoir  level  change  which  at  Canyon  Ferry  is  usually  on  the  order  of  a  few 
feet  or  less.  Annual  fluctuations  of  the  reservoir  level  at  Canyon  Ferry  are 
tjrpically  on  the  order  of  15  feet  or  so.  Thus  the  neglect  of  bank  storage 
changes  in  the  continuity  equation  is  not  expected  to  materially  affect  the 
equality. 

Groundwater  seepage  was  also  neglected  in  this  study.  The  seepage  term, 
RG,  represents  groundwater  flow  from  or  to  the  reservoir  site.  It  is  separate 
from  the  groundwater  flow  associated  with  the  temporary  on-site  bank  storage 
previously  discussed. 

The  net  groundwater  flow  is  almost  certainly  to  the  reservoir  since 
foundation  grouting  at  the  dam  severely  limits  groundwater  flow  from  the  site. 
The  primary  avenue  of  groundwater  inflow  is  believed  to  be  the  alluvial 
aquifer  beneath  and  adjacent  to  the  Missouri  River.  For  an  aquifer  width  of  2 
miles,  an  aquifer  transmissivity  of  120,000  gal/day/ft,  and  a  gradient  of  10 
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ft/mile  (assumed,  but  reasonable  values),  the  groundwater  inflow  would  be 
about  220  A-ft/mo.  This  inflow  is  negligible  compared  to  the  reservoir 
surface  inflow  RI  which  averages  171,000  A-ft  even  during  the  low  flow  month 
of  August.  Thus  the  neglect  of  reservoir  seepage,  RG,  will  have  little  effect 
on  the  equality  expressed  by  the  continuity  relationship. 

The  continuity  equation  was  further  simplified  by  the  neglect  of  pre- 
cipitation on  and  evaporation  from  the  Canyon  Ferry  Reservoir  surface  (RP  and 
RE  respectively).  This  neglect  disturbs  the  continuity  equality  somewhat  more 
than  the  previous  ones. 

Because  RP  and  RE  have  opposite  signs  (in  the  continuity  equation),  the 
neglect  of  one  compensates  in  part  for  the  neglect  of  the  other.  The  signifi- 
cance of  neglecting  RP  and  RE  varies  from  month  to  month,  but  is  generally 
greatest  in  August  when  evaporation  is  high  and  precipitation  is  low.  An 
evaporation  loss  of  5  inches  and  a  precipitation  gain  of  1  inch  (typical 
August  figures)  result  in  a  net  loss  of  approximately  12,000  A-ft  from  the 
35,000  acre  reservoir  surface.  This  is  7  percent  of  the  average  reservoir 
inflow  during  August  (171,000  A-ft).  The  neglect  of  evaporation  and  pre- 
cipitation during  other  months  is  relatively  less  significant  of  course, 
because  the  evaporation-precipitation  differential  is  generally  less  and 
because  the  reservoir  inflow  is  considerably  greater. 

On  the  basis  of  the  idealizations  discussed  above,  the  continuity 
equation  for  the  reservoir  takes  on  the  much  simpler  form 

RI  -  RG  =  RS   -  RS 
This  equation  provides  the  basis  for  operating  the  reservoir.   Recall  that  the 
known  (computed)  outflow  from  the  aggregate  irrigated  area  becomes  inflow  to 
Canyon  Ferry  (RI) .   Thus  for  a  specified  beginning-of-the-month  storage  RS., 
the  reservoir  outflow  RG  can  be  varied  in  such  a  wav  that  a  desired 
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end-of-the-month  storage  RS„  may  be  achieved.  The  monthly  outflow  and 
end-of-the-month  storage  are  generally  selected  to  meet  some  reservoir  opera- 
tional objective  or  objectives. 

RESERVOIR  OPERATION  PROCEDURES 

Canyon  Ferry  is  a  multipurpose  facility.  In  addition  to  its  power 
production  function,  it  provides  flood  control,  irrigation  water  storage,  and 
recreation.  It  also  serves  an  important  function  in  the  maintenance  of  flows 
for  downstream  power  producers  and  recreationalists.  All  of  these  functions 
(and  possibly  others)  should  be  considered  in  the  operation  of  the  reservoir. 
That  is  not  to  say  that  they  should  all  be  given  equal  consideration.  Some 
functions  are  more  important  than  others.  Some  functions  are  satisfied  by 
default.  Reservoir  recreational  objectives,  for  example,  are  reasonably  met 
within  the  reservoir  surface  level  range  dictated  by  flood  control  and  power 
generation  considerations. 

The  reservoir  operation  procedures  (rules)  used  in  the  simulations  of 
this  study  are  presented  and  discussed  below.  Operational  decisions,  based  on 
these  rules,  are  made  on  a  monthly  basis. 

1.  Reservoir  capacity  (2,050,900  A-ft  at  water  surface  elevation  3,800 
ft)  cannot  be  exceeded.  Water  is  released  (regardless  of  downstream 
channel  capacity)  to  meet  this  requirement. 

2.  Downstream  flow  is  limited  to  channel  capacity  (15,000  cfs)  whenever 
reservoir  storage  is  available  to  so  limit. 

3.  Reservoir  storage  in  excess  of  1,946,624  A-ft  (water  surface  ele- 
vation 3,797  ft)  is  reserved  for  flood  control.  This  storage  is 
used  only  when  failure  to  do  so  would  cause  the  downstream  channel 
capacity  (15,000  cfs)  to  be  exceeded. 

4.  Reservoir  releases  for  Helena  Valley  irrigation  and  pump  turbine 
discharge  are  maintained  at  historical  levels  regardless  of  the 
impact  on  reservoir  storage. 


26 


5.  Reservoir  releases  to  satisfy  the  minimum  downstream  flow  require- 
ment (2,500  cfs)  are  made  regardless  of  the  impact  on  reservoir 
storage. 

6.  When  not  limited  by  the  constraints  specified  above  (1-5),  reservoir 
releases  are  dictated  by  specified  month-end  storage  goals.  The 
amount  of  water  released  is  such  that  the  month-end  storage  goal 
will  be  met  (or  at  least  approximated) . 

7  Subject  to  the  constraints  above  (1-6),  and  the  performance  charac- 
teristics of  the  turbines,  turbine  discharge  is  regulated  in  a 
manner  designed  to  maximize  energy  production. 


Operation  rules  1-3  require  little  discussion.  Observance  of  rule  1 
assures  that  the  dam  will  not  be  overtopped.  Rule  2  requires  that  any  avail- 
able storage  must  be  used  to  reduce  downstream  flows  to  or  toward  15,000  cfs. 
Rule  3  requires  that  the  top  3  ft  of  reservoir  storage  be  reserved  exclusively 
for  flood  control.  This  storage  is  used  only  when  failure  to  do  so  would 
cause  flooding  (discharges  greater  than  15,000  cfs)  downstream. 

Rule  4  recognizes  the  priority  of  the  Helena  Valley  and  pump  turbine 
reservoir  releases.  Actually  these  releases  have  historically  been  quite 
limited  and  consequently  have  little  effect  on  reservoir  operation. 

A  minimum  downstream  flow  is  provided  by  Rule  5.  Thus  downstream  users 
(power  producers,  irrigators,  recreationalists,  etc.)  are  assured  a  water 
supply  during  periods  of  water  shortage.  Although  somewhat  arbitrary,  the 
minimum  flow  requirement  used  in  this  study  (2,500  cfs)  is  believed  to  be 
quite  reasonable. 

Rule  6  addresses  the  seasonal  differences  In  reservoir  operation  proce- 
dures that  are  mandated  by  seasonal  changes  in  reservoir  Inflow  and  user 
needs.  The  rule  is  applied  within  the  constraints  of  the  previous  5  rules. 
During  some  months,  reservoir  operation  is  totally  decerrained  by  rules  1-5. 
For  example,  downstream  flow  is  limited  to  channel  capacity  whenever  reservoir 
storage  Is  available  according  to  rule  2.   Furthermore,  if  the  only  available 


storage  is  flood  control  storage,  rule  3  requires  that  the  downstream  channel 
be  used  to  capacity.  Thus  releases  are  totally  constrained  and  rule  6  is  not 
considered.  However,  during  most  months,  operation  is  not  completely  con- 
strained by  rules  1-5,  and  it  is  during  these  months  that  rule  6  becomes 
operational. 

Rule  6  requires  that  the  month-end  reservoir  storage  be  brought  to  (or 
near)  a  specified  storage  goal  provided  rules  1-5  are  not  violated  in  so 
doing.   This  Is  accomplished  by  increasing  or  decreasing  reservoir  releases. 

The  month-end  storage  goals  are  based  on  historically  observed  month-end 
storage  volumes  at  Canyon  Ferry.  The  January,  February,  etc.  storage  goals 
are  essentially  the  averages  of  the  January,  February,  etc.  month-end  storage 
volumes.  Storage  goals  range  from  1,500,000  A-ft  (water  surface  elevation 
3,783  ft)  in  March  to  1,940,000  A-ft  (water  surface  elevation  3,797  ft)  in 
June. 

Rule  7  requires  that  the  available  releases,  as  determined  above  (1-6), 
be  used  to  maximize  energy  production.  In  other  words,  water  must  not  be 
spilled  at  the  expense  of  energy  production. 

It  should  be  noted  that  these  seven  reservoir  operation  rules  only 
approximate  actual  reservoir  operation  procedures.  They  do  not,  for  example, 
consider  the  recreational  use  of  the  reservoir.  However,  as  was  previously 
pointed  out,  this  use  is  reasonably  satisfied  within  the  constraints  imposed 
by  rules  1-7. 

The  seven  reservoir  operation  rules  outlined  and  discussed  above  ade- 
quately account  for  the  most  important  factors  upon  which  Canyon  Ferry  op- 
erations depend.  They  do  so  in  several  ways.  Most  reservoir  use  objectives 
such  as  flood  control,  energy  production,  etc.  are  of  course  directly  ad- 
dressed by  the  rules.   Other  objectives,  such  as  the  maintenance  of  reservoir 
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recreational  opportunities,  though  not  directly  addressed  by  the  rules,  are 
reasonably  well  met  within  the  constraints  imposed  by  the  rules.  Addition- 
ally, the  reservoir  storage  goals  of  rule  6  play  a  not  so  obvious  role  In 
meeting  operational  objectives.  These  goals  are  based  on  historical  storage 
records.  The  historical  storage  volumes  are  the  results  of  decisions  which 
were  presumably  aimed  at  meeting  the  combined  reservoir  use  objectives.  The 
goals  therefore  implicitly  reflect  all  of  the  combined  operational  objectives. 
It  follows  that  all  operational  objectives  are  at  least  indirectly  addressed 
when  these  goals  are  targeted  in  rule  6. 

Although  quite  reasonable  overall,  the  reservoir  operation  procedures  do 
have  certain  deficiencies  which  should  be  noted.  In  actual  reservoir  op- 
erations, snowmelt  runoff  forecasts  are  considered  in  the  determination  of 
spring  reservoir  releases.  These  forecasts  are  not  utilized  in  the  simu- 
lations of  this  study.  However,  the  snowmelt  management  objective  is  indi- 
rectly considered  in  that  the  historically  based  month-end  storage  goals 
reflect  that  objective. 

A  second  operational  deficiency  results  from  the  fact  that  the  rules  are 
applied  on  a  monthly  basis.  Actual  reservoir  operation  decisions  may  be  made 
on  a  daily  (or  even  hourly)  basis.  Although  otherwise  generally  adequate,  the 
monthly  operations  schedule  is  inadequate  for  flood  management.  A  flood  of 
several  days  duration  obviously  cannot  be  properly  managed  with  a  monthly 
operations  schedule.  The  inability  to  manage  floods  is  of  limited 
significance  however,  given  the  infrequency  of  occurrence  of  floods  that 
require  reservoir  management. 

One  other  operational  deficiency  results  from  the  rigidity  of  the  rules. 
Actual  reservoir  operation  procedure  are  more  flexible  and  sometimes  change 
over  time.   A  perusual  of  the  Canyon  Ferry  operation  records,  for  instance. 
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suggests  that  energy  production  has  become  an  increasingly  more  important 
operations  objective  over  the  years.  It  would  also  appear  that  the 
reservoir's  flood  control  function  has  changed  over  the  life  of  the  facility. 
Such  changes,  particularly  progressive  changes,  in  operational  objectives  are 
very  difficult  to  simulate  with  rigid  reservoir  operations  rules.  In- 
flexibility in  reservoir  operation  procedures  is  not  totally  disadvantageous 
however.  The  effects  of  changes  in  irrigation  activity,  for  example,  can  be 
better  assessed  when  viewed  in  the  framework  of  an  unchanging  reservoir 
operations  procedure. 

The  operational  deficiencies  described  in  the  previous  paragraphs,  though 
disadvantageous,  are  not  seriously  so.  Thus  the  overall  operation  of  the 
reservoir  can  be  satisfactorily  simulated  with  the  seven  specified  reservoir 
operation  rules. 

METHODOLOGY  SUMMARY 

Simulation  of  the  system  of  Figure  I  Is  accomplished  in  the  following 
manner.  The  hydrologlc  response  of  the  aggregate  irrigated  area  is  determined 
for  a  specified  sequence  of  irrigation  diversions.  The  continuity  equation 
and  the  ground  water  return  flow  equation  are  the  key  elements  in  this  re- 
sponse determination.  The  response  thus  determined,  the  irrigated  area 
outflow,  becomes  the  inflow  to  Canyon  Ferry  Reservoir.  That  inflow  is  routed 
through  the  reservoir  in  a  manner  that  satisfies  the  continuity  equation  and 
also  satisfies  (to  the  extent  possible)  the  reservoir  use  objectives. 

Application  of  the  methodology  to  the  Canyon  Ferry  system  is  described  in 
the  following  section. 
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CANYON  FERRY  SYSTEM 

Canyon  Ferry  Reservoir  is  located  on  the  Missouri  River  in  west  central 
Montana  (Figure  10).  The  reservoir  is  a  multipurpose  facility.  Among  its 
most  important  purposes  are  storage,  regulation  of  downstream  flows,  flood 
control,  energy  production,  and  recreation.  The  reservoir  has  a  surface  area 
of  35,000  acres  and  a  total  storage  volume  of  2,050,000  Acre-ft. 

The  Canyon  Ferry  catchment  includes  most  of  southwest  Montana.  It  is 
bounded  by  the  continental  divide  to  the  west  and  south  and  by  the  Gallatin 
and  Bridger  ranges  to  the  east.  The  major  streams  that  drain  the  area  are  the 
Gallatin,  Madison,  and  Jefferson  Rivers,  plus  the  two  major  Jefferson  River 
tributaries,  the  Big  Hole  and  the  Beaverhead  Rivers.  The  Missouri  River  is 
formed  by  the  confluence  of  the  Gallatin,  Madison,  and  Jefferson  Rivers. 
Southwestern  Montana  is  mountainous  and  the  streams  are  fed  primarily  by 
snowmelt  with  peak  discharges  generally  occurring  in  May  or  June. 

The  climate  in  southwestern  Montana  is  simiarid  with  annual  precipitation 
accumulations  on  the  order  of  15  inches  or  less  in  the  tillable  valleys. 
Where  available,  irrigation  waters  are  used  to  augment  precipitation  in  the 
production  of  hay  and  other  crops. 

The  features  of  the  Canyon  Ferry  system  that  are  relevant  to  the  method- 
ology of  this  study  are  discussed  under  separate  headings  below. 
AGRICULTURAL  ACTIVITIES 

Although  irrigation  is  practiced  at  many  locations,  probably  90  percent 
of  it  is  concentrated  at  the  locations  identified  by  the  cross-hatched  areas 
of  Figure  10.  Irrigation  activities  are  in  general  quite  similar  at  the 
various  locations.  Irrigation  waters  are  diverted  from  streams  and  applied  to 
lands  chat  are  in  close  proximity  to  the  streams  (bottom  lands,  benches,  and 
alluvial  fans).   Though  occupying  different  stream  basins,  the  irrigated  lands 
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Figure  10.   Ceinyon  Ferry  study  area. 
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are  of  similar  geologic  origin  and  depositional  character,  and  thus 
hydrologlcaliy  similar.  Irrigation  methods  ar«  also  quite  similar  throughout 
the  drainage.  Flood  irrigation  using  contour  ditches  to  deliver  the  diverted 
surface  water  to  the  fields  is  most  common.  Sprinkler  Irrigation,  though  less 
common,  has  become  increasingly  more  Important  in  recent  years. 

As  discussed  at  length  in  the  methodology  section,  these  irrigated  areas 
(hydrologlcaliy  similar  and  similarly  managed)  are  treated  in  aggregate  in 
this  study. 

Irrigation  Diversions  and  Crop  Consumptive  Use 

Irrigation  waters  are  diverted  from  streams  at  many  locations  within  the 
basin.  Unfortunately,  the  diversion  flow  rates  are  regularly  measured  at  only 
a  very  few  locations.  Diversion  totals  (and  consumptive  use  totals)  for  the 
basin  must  therefore  be  estimated.  The  Water  Conservation  and  Salvage  Report 
for  Montana,  a  Soil  Conservation  Service  publication  (1),  provides  such 
estimates.  The  estimates  are  based  on  available  diversion  measurements,  data 
on  crop  types,  acreages,  and  consumptive  needs,  and  on  firsthand  knowledge  of 
the  areas  and  the  irrigation  operations. 

The  report  estimates  annual  irrigation  diversions,  crop  consumption 
needs,  etc.  on  county  by  county  basis.  The  monthly  basin-wide  figures 
required  for  this  study  were  obtained  from  the  annual  county  figures  in  the 
following  manner.  Annual  basin-wide  totals  are  determined  from  the  county 
figures  and  from  the  hydrologic  maps  in  the  publication  titled  "An  Atlas  of 
Water  Resources  in  Montana  by  Hydrologic  Basins"  (A) .  The  maps  show  the  basin 
boundaries,  the  county  boundaries,  and  the  Irrigated  areas  within  the  counties 
(and  basin).   They  thus  provide  a  basis  for  determining  what  fraction  of  each 
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county's  irrigation  activity  takes  place  within  the  basin  and  correspondingly 
what  fraction  of  that  county's  diversions,  etc.  are  within  the  basin.  The 
annual  basin  totals  (the  summed  county  fractions)  are  then  allocated  on  a 
monthly  basis  over  the  irrigation  season.  The  allocation  is  based  on  histor- 
ical diversion  records  from  the  region.  The  presumption  is  that  these  records 
are  indicative  of  diversion  timing  throughout  the  basin. 

The  annual  basin  wide  totals  for  irrigation  diversion  and  for  crop 
consumptive  use  were  found  to  be  3,503,400  Acre-ft  and  651,300  Acre-ft 
respectively.  Diversion  timing  (and  consumptive  use)  was  assumed  to  be  as 
follows:  25%  of  the  annual  total  diverted  in  each  of  the  three  months  of 
June,  July,  and  August;  lOZ  of  the  annual  total  in  each  of  the  months  of  May 
and  September;  and  2.5%  of  the  annual  total  in  each  of  the  months  of  April  and 
October.  The  amount  and  timing  of  the  diversions  (and  consumptive  use)  were 
assumed  to  be  the  same  from  year  to  year.  Additional  details  about  the 
determination  of  the  basin-wide  irrigation  figures  can  be  found  in  Appendix  A. 

Irrigation  Efficiency 

Irrigation  efficiency  Is  herein  defined  as  the  ratio  of  the  crop 
consumptive  use  to  the  irrigation  diversion.  Based  on  the  annual  totals  of 
the  previous  paragraph,  it  follows  that  the  irrigation  efficiency  of  the 
aggregate  area  is  19%  (651,300  t  3,503,400).  The  irrigation  efficiency  is 
assumed  to  be  the  same  (19%)  for  each  month  of  the  irrigation  season. 

A  very  large  portion  of  the  irrigation  diversion  (81%)  is  not  beneficial- 
ly used  in  crop  production.  The  disposal  of  this  water  is  discussed  in  the 
following  sections. 
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Nonbenef iclal  Consumptive  Use 

The  Water  Conservation  and  Salvage  Report  for  Montana  (1)  states  that  in 
addition  to  beneficial  crop  consumptive  use  within  the  basin,  there  is  also 
nonbenef icial  consumptive  use.  This  nonbenef icial  use  (by  phreatophyes  that 
border  ditches,  etc.)  is  estimated  to  be  25Z  of  the  beneficial  use.  This 
water,  along  with  the  beneficially  consumed  water,  is  of  course  lost  to  the 
Canyon  Ferry  hydrologic  system. 
Surface  Return  Flow 

A  portion  of  the  unconsumed  irrigation  diversion  returns  to  the  stream 
via  the  surface.  During  the  months  of  April,  May,  and  June,  it  is  assumed 
that  25Z  of  each  month's  diversion  returns  this  way.  During  July,  August, 
September,  and  October,  15%  of  each  month's  diversion  is  assumed  to  return  in 
this  way.  Although  somewhat  arbitrary,  these  percentages  are  believed  to  be 
reasonable,  and  to  reflect  to  some  degree  the  more  relaxed  management  that 
prevails  during  the  water  rich  months  of  April,  May  and  June. 

It  is  assumed  that  these  surface  flows  return  to  the  stream  during  the 
same  month  in  which  they  are  diverted.  This  is  quite  reasonable  in  view  of 
the  close  proximity  of  the  irrigated  areas  to  the  streams. 

Groundwater  Recharge 

The  remainder  of  the  irrigation  diversion  (the  irrigation  diversion  minus 
the  consumed  water  and  the  surface  return)  recharges  the  ground  water  aquifer. 
It  does  so  via  ditch  seepage  losses  and  deep  percolation  through  the  root 
zone.  It  is  assumed  that  the  groundwater  system  is  recharged  during  the  same 
month  in  which  the  water  is  diverted  and  applied  to  the  fields.   This 
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assumption  does  not  imply  that  the  downward  moving  water  necessarily  traverses 
the  partially  saturated  zone  in  a  very  short  period  of  time,  but  rather  it 
implies  that  the  soil  moisture  content  in  the  partially  saturated  zone  is 
constant  (near  field  capacity)  and  consequently  recharge  (excess  over  crop 
consumptive  need)  to  that  zone  from  above  must  equal  discharge  from  that  zone 
to  the  ground  water  system  below. 

GROUNDWATER 

Ground  water  return  flows  move  from  the  recharge  areas  (irrigated  areas) 
through  the  aquifer  to  the  surface  drainage  system.  The  time  delay  associated 
with  that  movement  is  dependent  upon  the  aquifer  properties  and  upon  the 
location  and  configuration  of  the  drains. 

Aquifer  Properties 

In  the  early  1950's,  the  United  States  Geological  Survey  (USGS)  carried 
out  a  comprehensive  study  of  the  groundwater  resources  of  the  Gallatin  Valley 
(5).  As  part  of  the  study,  aquifer  tests  were  conducted  to  determine  the 
properties  of  the  alluvial  deposits  underlying  the  irrigated  areas  of  that 
valley.  These  test  results  provided  the  basis  for  quantifying  the  aggregate 
aquifer  properties.  Specifically,  it  was  assumed  that  the  aggregate  aquifer 
properties  were  the  same  as  those  of  the  Gallatin  Valley  aquifer.  Although 
somewhat  of  a  necessity  (no  other  substantive  aquifer  tests  have  ever  been 
conducted  in  the  areas  of  concern) ,  the  assumption  is  never-the-less  believed 
to  be  quite  reasonable.  The  irrigated  lands  and  underlying  aquifers  are 
hydrologically  similar.   Though  in  different  basins,  they  typically  occupy 
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similar  settings  close  to  streams.  The  aquifers  are  also  of  similar  geologic 
origin  and  depositional  character. 

In  the  course  of  the  U.S.  Geological  Survey  study,  aquifer  tests  were 
conducted  at  numerous  locations  across  the  irrigated  areas  of  the  Gallatin 
Valley.  Test  results  indicate  fairly  large  variations  in  aquifer  properties 
with  little  or  no  spatial  pattern  in  those  property  variations.  This  is  not 
surprising  in  view  of  the  depositional  character  of  the  aquifer.  Sediments 
have  been  carried  from  the  nearby  mountains  by  numerous  streams.  Over  geolog- 
ic time,  these  streams  have  changed  continuously  in  terms  of  discharge, 
sediment  load,  channel  configuration,  and  channel  location.  The  result  is  an 
alluvial  aquifer  that  exhibits  a  tremendous  variation  in  sediment  size, 
ranging  from  large  boulders  to  silt  and  clay  size.  Though  generally  well 
mixed,  there  are  deposits  of  well  sorted  materials  such  as  sand  or  gravel 
stringers  and  clay  lenses.  These  well  sorted  deposits  are  typically  limited 
in  extent  and  more-or-less  randomly  distributed. 

The  alluvial  deposits  described  above  are  of  Quaternary  age.  They  are 
underlain  by  much  less  permeable  deposits  of  Tertiary  age.  The  alluvial 
deposits,  the  concern  of  this  study,  are  variable  in  thickness.  The  USGS 
study  found  alluvium  thickness  to  vary  from  only  a  few  feet  to  as  much  as  400 
feet  with  thicknesses  on  the  order  of  50  to  120  feet  being  most  common. 

The  alluvial  Quaternary  aquifer  described  in  the  previous  paragraphs  is 
obviously  not  homogeneous  -  at  least  not  on  a  local  scale.  This  is  borne  out 
by  the  variability  in  the  USGS  aquifer  property  values.  However,  the  absence 
of  any  pattern  in  the  spatial  variation  in  aquifer  properties  suggests  that 
the  aquifer  might  be  considered  homogeneous  on  a  large  scale.  In  other  words, 
the  randomly  heterogeneous  aquifer  with  its  varied  properties  may  be  rea- 
sonably approximated  by  a  homogeneous  aquifer  with  properties  that  are  the 
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average  of  the  hecerogeneous  aquifer  properties.  This  sort  of  Idealization  Is 
appropriate  when  Investigating  the  overall  behavior  of  the  Irrigated  area 
ground  water  system.  It  would  not  be  appropriate  for  Investigating  localized 
ground  water  flows  within  the  Irrigated  area.  Because  the  concern  of  this 
study  Is  overall  aquifer  behavior,  the  aquifer  Is  considered  to  be  homogeneous 
vlth  averaged  aquifer  properties.  The  assumption  of  homogenlety  employed  In 
the  development  of  the  return  flow  equation  (methodology  section)  Is  thus 
justified. 

A  condition  of  aquifer  Isotropy  was  also  assumed  in  the  development  of 
Che  return  flow  relationship.  Unfortunately,  no  aquifer  tests  were  conducted 
to  assess  the  Isotropy  of  the  Gallatin  Valley  aquifer.  The  mode  of  sediment 
deposition  la  the  alluvial  aquifer  however,  does  suggest  certain  things  about 
the  Isotropy  of  the  aquifer.  When  equldlmenslonal  sediments  are  deposited  by 
flowing  water,  they  have  no  preferred  deposltional  orientation.  However,  when 
nonequldlmenslonal  particles  such  as  those  most  commonly  found  In  Intermoun- 
taln  alluvial  basins  are  deposited,  they  tend  to  be  positioned  with  the 
largest  dimensions  parallel  to  the  deposltional  surface  (Figure  II). 


Figure    11.     Deposltional    orientation    of   alluvial    sediments. 
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This  preferred  orientation  causes  the  deposits  to  be  anisotropic  with  the 
hydraulic  conductivity  in  the  vertical  z  direction  typically  significantly 
less  than  the  equal  horizontal  x  and  y  direction  conductivities.  The 
anisotropy  must  be  taken  into  account  when  the  ground  water  flow  velocity  has 
a  significant  vertical  component.  However,  when  the  flow  is  strictly  horizon- 
tal, the  anisotropy  is  not  a  factor.  The  ground  water  flow  in  the  irrigated 
areas  is  predominantly  horizontal,  and  consequently  the  assumption  of  aquifer 
isotropy  in  the  return  flow  equation  development  is  not  unreasonable.  The 
horizontal  nature  of  the  irrigated  area  groundwater  flows  is  discussed  in  more 
detail  in  the  next  section  on  drain  configuration. 

USGS  tests  conducted  within  the  irrigated  areas  of  the  Gallatin  Valley 
yielded  values  for  aquifer  transmissivity  T  (aquifer  conductivity  times 
aquifer  saturated  thickness)  ranging  from  4,500  to  670,000  gal/day/ft. 
Transmissivities  were  generally  much  less  variable  than  these  extremes  might 
suggest,  however.  The  average  transmissivity  on  the  alluvial  plain  (24  test 
sites)  was  212,000  gal/day/ft.  The  average  transmissivity  on  the  alluvial  fan 
deposits  (6  test  sites)  was  41,000  gal/day/ft.  The  average  aquifer 
transmissivity,  based  on  an  approximate  area  weighting,  was  assumed  to  be 
120,000  gal/day/ft.  The  individual  well  transmissivities  upon  which  the 
average  transmissivity  is  based  are  presented  in  tabular  form  in  Appendix  B. 

The  USGS  study  unfortunately  provides  only  limited  information  about  the 
storage  coefficient  S  of  the  Gallatin  Valley  aquifer.  The  aquifer  tests 
conducted  in  the  study  were  either  single  well  tests,  which  yield  no 
information  about  the  storage  coefficient,  or  short  duration  multiple  well 
tests.  Although  the  storage  coefficient  can  be  determined  from  a  multiple 
well  test,  the  coefficient  value  tends  to  be  deceptively  small  when  the  test 
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is  of  short  duration.  This  is  primarily  because  the  analysis  upon  which  the 
coefficient  determination  is  based  assumes  an  instantaneous  and  total  release 
of  available  water  from  the  pore  space  as  the  water  table  falls  during 
testing.  In  actuality,  the  available  water  is  released  at  an  exponentially 
decreasing  rate  over  time.  As  a  consequence  of  this  gradual  release,  the 
computed  value  of  the  storage  coefficient  tends  to  increase  and  approach  a 
constant  value  as  the  duration  of  the  aquifer  test  increases.  Because  this 
study  considers  ground  water  level  fluctuations  on  a  monthly  basis,  it  is 
concerned  with  the  longer  term  (1  month)  release  of  water  from  the  dewatered 
portion  of  the  aquifer  and  the  storage  coefficient  that  characterizes  that 
release. 

The  USGS  estimated  that  the  average  storage  coefficient  for  the  Gallatin 
Valley  alluvium  would  be  about  0.05  after  12  hours  of  pump  testing  and  much 
larger  after  several  weeks  of  pumping.  On  the  basis  of  this  information,  and 
on  the  basis  published  storage  coefficient  estimates  for  similar  alluvial 
deposits,  a  coefficient  of  0.20  was  assumed  appropriate  for  the  aquifer. 

These  values  of  the  storage  coefficient  S  and  the  transmissivity  T  are 
used  in  the  determination  of  the  groundwater  return  flows  from  the  Canyon 
Ferry  aggregate  irrigated  area.  Recall  from  the  methodology  section  that  the 
return  flow  determination  was  based  on  the  return  flow  factors,  r(l,j),  and 
that  these  factors  are  functions  of  — rr,  where  a  =  T/S,  t  is  time,  and  L  is 
drain  spacing  (discussed  in  next  section),  A  tabulation  of  the  Canyon  Ferry 
system  return  flow  factors  is  presented  in  Appendix  C. 
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Drala  Configuration 

The  groundwater  return  flow  determination  procedures  developed  in  the 
methodology  section  were  for  an  idealized  surface  drain  configuration.  The 
procedures  are  appropriate  for  a  system  of  parallel  drains  spaced  a  distance  L 
apart  or  for  a  single  drain  centered  in  a  valley  of  width  L  (Figure  8) .  These 
Idealizations  may  be  only  approximately  satisfied  in  real  groundwater  drainage 
systems. 

Maps  prepared  by  the  Montana  State  Engineer's  Office  (6,7)  show  that 
irrigated  lands  in  southwestern  Montana  are  traversed  by  numerous  effluent 
streams  and  drains.  Although  the  configurations  of  these  drains  (natural  and 
man-made)  are  varied,  they  can  generally  be  roughly  approximated  by  the 
Idealized  configurations  described  above.  The  configuration  of  a  single  drain 
centered  in  a  valley  of  width  L  is  approximated  in  some  of  the  narrower 
irrigated  valleys  in  the  Canyon  Ferry  drainage  (Figure  12). 

Stream 


A' 


Section  A-A' 


Figure  12.   Valley-centered  drain  (scream) 
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lo  the  broader  iQCennouncalo  alluvial  basins  such  as  Che  Gallacla  Valley,  Che 
surface  drainage  network  In  the  Irrigated  areas  tends  to  be  dendritic.  These 
dendritic  networks  closely  resemble  a  system  of  parallel  lateral  drains 
connected  by  mains  in  a  herringbone  pattern  (Figure  13). 


Scream 


Mains 


Figure  13.  Dendritic  drain  pattern. 
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Section  B-B' 


Thus  the  areal  configurations  of  the  existing  drains  in  the  Canyon  Ferry 
drainage  are  reasonably  approximated  by  combinations  of  the  idealized  drain 
systems  of  the  methodology  section.  The  distance  between  "parallel"  drains, 
as  determined  from  the  State  Engineer's  Office   publications   (6,7)   is 
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estimated  to  average  about  1.5  to  2  miles.   A  drain  spacing  L  of  2  miles  was 
assumed  in  this  study. 

In  addition  to  the  Idealized  areal  drain  configuration,  certain  other 
idealizations  were  employed  in  the  development  of  the  groundwater  return  flow 
computation  procedures  in  the  methodology  section.  Specif Ically, it  was 
assumed  that  the  drainable  depth  h  was  much  less  than  the  saturated  thickness 
D  (Figure  4),  that  the  groundwater  flow  was  horizontal  and  one-dimensional  (in 
the  X  direction  of  Figure  4) ,  and  that  the  Dupuit-Forchheimer  Idealization  was 
appropriate  for  characterizing  the  groundwater  flow  velocity.  The  appropri- 
ateness of  these  three  assumptions  is  dependent  upon  the  saturated  thickness 
of  the  aquifer,  the  uniformity  of  that  thickness,  and  the  vertical  position  of 
the  drain  in  the  saturated  thickness.  As  with  the  previously  discussed  areal 
drain  configuration  idealization,  these  idealizations  may  be  only  approximated 
in  the  real  groundwater  system.  With  reference  to  the  typical  aquifer 
sections  of  Figures  12  and  13  (sections  A-A'  and  B-B'  respectively),  it  can  be 
seen  that  the  drainable  depth  is  much  less  than  the  aquifer's  saturated 
thickness.  Recall  that  the  aquifer  thickness,  although  somewhat  variable,  was 
generally  in  excess  of  50  feet.  The  depth  of  the  drains  is  generally  on  the 
order  of  10  feet  or  less,  resulting  in  a  drainable  depth  of  less  than  10  feet 
and  a  saturated  thickness  in  excess  of  40  feet.  With  reference  again  to  the 
cross  sections  of  Figures  12  and  13,  it  can  also  be  seen  that,  though  not 
strictly  horizontal,  the  groundwater  flow  is  predominately  horizontal  and 
one-dimensional  in  a  direction  toward  the  drain.  The  appropriateness  of  the 
Dupuit-Forchheimer  assumption  is  also  apparent  from  the  cross  sections  in  that 
the  slope  of  the  water  table  is  small  and  the  aquifer's  saturated  thickness  is 
large  as  required  for  application  of  the  assumption. 
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Thus  the  idealizations  employed  in  the  development  of  the  groundwater 
return  flow  computation  procedures  in  the  methodology  section  are  reasonably 
satisfied  in  the  real  groundwater  system.  It  follows  that  those  computational 
procedures  can  reasonably  be  used  to  determine  groundwater  return  flows  from 
the  aggregate  irrigated  area  in  the  Canyon  Ferry  drainage. 

SURFACE  WATER 

A  greatly  simplified  surface  water  system  results  from  the  previously 
discussed  aggregation  of  Irrigated  areas  (Figure  1).  As  suggested  by  the 
figure,  streamflows  of  concern  are  those  flowing  into  the  irrigated  area, 
those  flowing  out  of  the  irrigated  area  (and  into  the  reservoir)  ,  and  those 
flowing  out  of  the  reservoir.  Also  of  concern  of  course, is  the  temporary 
storage  afforded  the  streamflow  by  the  Canyon  Ferry  Reservoir. 

Streamflow 

Recall  that  this  study  considers  total  inflow  to  and  total  outflow  from 
the  aggregate  irrigated  area  (Figure  1).  It  does  not  consider  individual 
stream  inflows  or  outflows.  Recall  also  that  the  total  inflow,  the  "natural" 
Inflow,  is  computed  from  the  known  total  outflow  by  means  of  the  continuity 
and  return  flow  relationships.  The  "natural"  inflows  thus  computed  then 
become  the  surface  water  data  base  for  subsequent  analyses  as  the  inflows  are 
subjected  to  various  assumed  irrigation  use  demands  and  the  basin  response  is 
simulated  (irrigated  area  outflows,  etc.  are  simulated).  Thus  the  only 
streamflow  records  required  for  the  desired  hydrologic  simulations  are  those 
for  total  outflow  from  the  aggregate  irrigated  area.   Total  outflow  is 
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essentially  the  same  as  the  Inflow  to  Canyon  Ferry  Reservoir  (see  discussion 
under  "Aggregation  of  Irrigated  Areas"  section),  and  records  of  Canyon  Ferry 
inflow,  as  determined  by  the  Bureau  of  Reclamation,  are  available  in  the 
Bureau's  Canyon  Ferry  Reservoir  operation  records  (Appendix  D)  .  Reservoir 
inflow  records  from  1955  (the  initial  filling  of  the  reservoir)  through  1984 
were  utilized  in  this  study. 

Canyon  Ferry  Reservoir 

Outflows  from  the  aggregate  irrigated  area  enter  Canyon  Ferry  Reservoir, 
are  temporarily  stored  in  the  reservoir,  and  are  eventually  released  from  the 
reservoir.  Operation  of  the  reservoir  (release  scheduling,  etc.)  depends  on 
many  factors.  The  numerous  use  objectives  of  the  multipurpose  facility  must 
of  course  be  considered.  Physical  system  constraints,  such  as  reservoir 
capacity  and  turbine  capacity,  must  also  be  considered  along  with  operational 
constraints  such  as  minimum  release  requirements.  These  and  other  factors 
that  relate  to  the  Canyon  Ferry  Reservoir  operations  are  discussed  in  the 
following  paragraphs. 

Various  in-house  documents  from  the  Bureau  of  Reclamation  describe  the 
Canyon  Ferry  facility,  its  operational  objectives,  and  its  historical  op- 
erations. Information  from  these  documents  provided  the  framework  upon  which 
the  simulation  of  reservoir  operations  was  based.  Of  particular  importance 
was  the  Canyon  Ferry  Reservoir  operations  record  (Appendix  D) .  This  document 
lists  reservoir  inflows,  outflows,  storage  volumes,  etc.  on  a  monthly  basis 
over  the  life  of  the  reservoir. 
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Water  may  be  released  from  the  Canyon  Ferry  Reservoir  via  the  river  and 
spillway  outlets,  through  the  power  turbines,  through  the  pump  turbine,  or  it 
may  be  pumped  to  the  Helena  Valley.  In  the  simulations  of  this  study,  re- 
leases through  the  river  and  spillway  outlets  and  through  the  power  turbines 
are  allowed  to  vary.  They  are  adjusted  (via  reservoir  operations)  in  response 
to  changes  in  reservoir  inflow.  Pumpage  to  Helena  Valley  and  pump  turbine 
discharge,  however,  are  not  allowed  to  vary.  These  releases  are  maintained  at 
the  historical  levels  (timing  and  amount)  indicated  in  the  Canyon  Ferry 
Reservoir  operations  record. 

As  discussed  in  detail  in  the  methodology  section  (under  the  "Reservoir 
Operation  Procedures"  heading) ,  monthly  reservoir  releases  are  determined  in 
part  by  month-end  storage  goals.  Within  the  limits  imposed  by  other  con- 
straints (minimum  and  maximum  release  requirements,  etc.),  releases  are 
Increased  or  decreased  so  that  the  month-end  storage  objectives  can  be  at- 
tained, or  at  least  approached.  The  month-end  goals,  which  are  intended  to 
reflect  historical  reservoir  operations  objectives,  are  obtained  essentially 
by  averaging  the  month-end  storage  volumes  recorded  in  the  Canyon  Ferry 
Reservoir  operations  record.  The  month-end  storage  goals  are  listed  in 
Appendix  E. 

In  addition  to  the  direct  use  of  the  reservoir  operations  records  de- 
scribed above,  the  records  were  also  used  as  a  basis  for  assessing  the  reason- 
ableness of  the  simulated  results.  This  was  accomplished  by  comparing  the 
various  simulated  results  (monthly  reservoir  releases,  storages,  etc.)  to 
their  historical  counterparts. 

Bureau  of  Reclamation  documents  relating  to  turbine  operating 
characteristics,  and  reservoir  storage  allocations  provide  additional 
information  that  is  utilized  in  the  simulation  of  reservoir  operations. 
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Simulated  power  turbine  performance  is  constrained  by  the  turbine  operating 
characteristics  as  defined  by  the  Canyon  Ferry  power  plant's  characteristic 
curves  (Appendix  F) .  Simulated  reservoir  operations  require  that  storage  be 
reserved  for  flood  control  as  specified  in  the  Bureau's  Canyon  Ferry  Reservoir 
allocations  (Appendix  F)  .  Total  reservoir  storage,  an  obvious  consideration 
in  the  simulation  of  reservoir  operations,  is  also  identified  in  this  latter 
document . 

A  Bureau  document  titled  "General  Operating  Principles  at  Canyon  Ferry" 
(Appendix  G)  was  extensively  used  in  the  formulation  of  rules  for  the  simu- 
lation of  reservoir  operations.  The  operating  principles  identify  various 
operational  objectives,  priorities,  and  constraints.  The  simulation  rules, 
described  in  detail  in  the  methodology  section,  adhere  to  these  principles  to 
the  extent  possible.  Deviations  from  these  principles  are  discussed  in  the 
methodology  section. 

Reservoir  releases  to  the  Missouri  River  are  limited  to  channel  capacity 
whenever  reservoir  storage  space  is  available  to  so  limit.  The  channel 
capacity  is  15,000  cfs  according  to  the  Bureau's  document  on  reservoir  operat- 
ing principles.  In  addition,  a  minimum  reservoir  release  of  2,500  cfs  is 
maintained  at  all  times  to  satisfy  downstream  users  (power  producers, 
recreationalists,  irrigators,  etc.).  Though  somewhat  arbitrary,  this  minimum 
flow  requirement  was  suggested  as  a  reasonable  limit  by  Mr.  Gordon  Aycock  of 
the  Bureau  of  Reclamation. 

SUMMARY  -  CANYON  FERRY  SYSTEM 

As  is  apparent  from  the  discussions  of  this  section,  it  has  been  neces- 
sary  to   use   numerous   simplifying   assumptions   and   idealizations   in 
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characterizing  the  Canyon  Ferry  system.  This  is  not  surprising  in  view  of  the 
size  and  diversity  of  the  system,  and  also  the  general  scarcity  of  data 
characterizing  the  system.  The  simplifications,  while  sacrificing  something 
in  detail  and  rigor,  are  essential  in  that  they  make  possible  a  reasonable 
simulation  of  the  total  basin's  hydrologic  behavior. 
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SIMULATIONS  AND  RESULTS 
Recall  that  the  primary  objective  of  this  study  is  to  assess  the 
hydrologic  impacts  that  result  from  changes  in  irrigation  activity  in  the 
Canyon  Ferry  drainage.  Toward  this  end,  numerous  simulations  of  the  basin's 
hydrologic  behavior  were  made.  The  first  simulation  yielded  a  record  of 
monthly  irrigated  area  inflow.  Subsequent  simulations  utilized  this  record  in 
assessing  the  Impacts  of  changes  in  irrigation  activity.  The  various  simu- 
lations are  discussed  on  the  following  pages. 

SIMULATION  OF  IRRIGATED  AREA  INFLOW 

As  described  in  detail  in  the  methodology  section,  irrigated  area  inflows 
are  determined  by  solving  the  continuity  equation  and  the  groundwater  return 
flow  equation  for  the  aggregate  irrigated  area.  The  "natural"  inflows  thus 
obtained  then  become  the  surface  water  data  base  for  subsequent  analyses  in 
which  these  inflows  are  subjected  to  various  assumed  irrigation  use  demands. 

Solution  of  the  continuity  and  return  flow  equations  for  irrigated  area 
inflow  requires  that  irrigated  area  outflows,  irrigation  diversions,  and 
irrigation  efficiencies  be  specified  along  with  drainage  system  parameters 
(aquifer  transmissivity ,  aquifer  storage  coefficient,  and  drain  spacing).  The 
variables  to  be  specified  are  briefly  discussed  in  the  following  paragraphs, 
and,  with  the  exception  of  irrigated  area  outflow,  are  tabulated  in  Table  1. 
A  tabulation  of  irrigated  area  outflows  (Canyon  Ferry  inflows)  may  be  found  in 
Appendix  D.  For  a  more  detailed  discussion  of  the  variables  and  procedures, 
refer  to  the  previous  sections  on  the  Canyon  Ferry  system  and  on  methodology. 

Irrigated  area  outflows  are  assumed  to  be  equal  to  Canyon  Ferry  inflows. 
Thus  "historical"  records  of  monthly  irrigated  area  outflow  may  be  determined 
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from  reservoir  inflows  as  recorded  in  the  Canyon  Ferry  operations  record 
(Appendix  D) . 

A  monthly  consumptive  use  schedule  (Table  1)  was  estimated  from  "The 
Water  Conservation  and  Salvage  Report  for  Montana"  (1). 

TABLE  1 
Simulation  Input  Data 
Crop  Consumptive  Use  Schedule 
Month Consumptive  Use  (Acre-ft) 


Oct 

Nov 

Dec 

Jan 

Feb 

Mar 

Apr 

May 

June 

July 

Aug 

Sept 


16,282.0 

0.0 

0.0 

0.0 

0.0 

0.0 

16,282.0 

65,130.0 

162,825.0 

162,825.0 

162,825.0 

65,130.0 


Irrigation  Efficiency  »  19% 

Aquifer  Transmissivity  -  120,000  gal/day/ft 

Aquifer  Storage  Coefficient  =  0.20 

Drain  Spacing  -  2  miles 


50 


This  report  also  yielded  an  irrigation  efficiency  (consumptive  use/irrigation 
diversion)  estimate  of  19  percent  (See  also  Appendix  A).  Both  the  consumptive 
use  schedule  and  the  irrigation  efficiency  were  assumed  to  be  the  same  from 
year  to  year  throughout  the  simulation  period.  These  assumptions  are  only 
approximately  true,  of  course.  They  were  made  out  of  necessity,  since  there 
are  no  reliable  records  that  tabulate  year  to  year  changes  in  irrigation 
activity.  Irrigation  diversions  were  computed  from  the  consumptive  use 
figures  and  the  known  irrigation  efficiency. 

The  return  flow  aquifer's  transmissivity  T  and  storage  coefficient  S  were 
estimated  to  be  120,000  gal/day/ft  and  0.20  respectively.  Spacing  between 
"parallel"  drains  L,  was  estimated  to  be  2  miles.  Groundwater  return  flow 
factors,  r(l,j),  are  a  function  of  these  variables  (T,  S,  and  L)  along  with 
time  t.  Return  flow  factors  were  determined  for  monthly  time  intervals  and 
are  presented  in  Appendix  C.  Recall  that  r(l,J)  is  the  percentage  of  a  unit 
groundwater  recharge  during  month  1  that  returns  to  the  surface  water  system 
(drains)  in  month  j.  The  factors  thus  provide  a  way  to  determine  the 
groundwater  return  flow  during  any  month  from  a  known  groundwater  recharge 
during  a  specified  month. 

Groundwater  recharge  during  a  particular  month  is  dependent  upon  the 
irrigation  diversion  during  that  month.  A  portion  of  the  diversion  is 
consumed,  a  portion  returns  via  the  surface  to  the  stream,  and  the  remainder 
recharges  the  groundwater  system. 

The  "Water  Conservation  and  Salvage  Report  for  Montana" (1)  estimates  that 
in  addition  to  the  crop  consumptive  use,  there  is  an  additional  consumptive 
use  (phreatophytes,  etc.)  that  is  equal  to  25  percent  of  the  crop  consumptive 
use.  Monthly  consumptive  use  was  therefore  assumed  to  be  123  percent  of  the 
consumptive  use  values  in  Table  1. 
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Surface  return  flows  from  irrigation  applications  typically  vary  through 
the  irrigation  season.  In  this  study  it  was  assumed  that  25  percent  of  the 
monthly  diversion  returns  to  the  stream  via  surface  runoff  in  the  water-rich 
months  of  April,  May,  and  June,  while  only  15  percent  returns  by  this  route  in 
the  water-short  months  of  July,  August,  September  and  October.  These  surface 
flows  return  during  the  same  month  that  they  are  diverted.  Irrigation  waters 
that  are  neither  consumed  nor  return  via  the  surface  constitute  groundwater 
recharge. 

In  summary  then,  from  the  consumptive  use  figures  of  Table  1,  and  the 
current  irrigation  efficiency  of  19  percent,  monthly  irrigation  diversions  can 
be  computed.  Monthly  surface  returns  are  taken  as  a  percentage  (specified 
above)  of  these  diversions.  Monthly  groundwater  recharge  is  found  by  sub- 
tracting the  monthly  surface  return  and  monthly  consumptive  use  from  the 
monthly  irrigation  diversion.  Monthly  groundwater  return  flow  is  in  turn 
determined  from  the  monthly  groundwater  recharge  via  the  return  flow  factors. 
And  finally,  irrigated  area  inflow  is  determined  from  the  continuity  equation, 

(I  +  TI  +  P  +  SRN  +  GN)  +  (SRI  +  GI)  -  D  -  0  »  0  (1) 

As  previously  discussed,  the  various  aggregate  area  inflows  (left  parentheses) 
need  not,  and  in  general  cannot  be  separately  identified.  Thus  the  continuity 
equation  can  be  written  as 

Inflow  +  (SRI  +  GI)  -  D  -  0  »  0 
where  inflow  is  total  irrigated  area  inflow.   The  equation  may  be  solved  for 
monthly  inflow  since  monthly  values  of  surface  runoff  from  irrigation, 
groundwater  return  from  irrigation,  diversion,  and  outflow  (SRI,  GI,  D,  and  0 
respectively)  are  known. 

Monthly  values  of  irrigated  area  inflow  were  determined  in  this  manner 
for  the  period  1955  -  1984.   These  "natural"  inflows  are  tabulated  in  Appendix 
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H,  and  a  typical  interval  from  that  tabulation  is  shown  graphically  as  a 
hydrograph  in  Figure  14.  The  irrigated  area  outflow  hydrograph  (Canyon  Ferry 
inflow  hydrograph)  from  which  the  "natural"  inflow  hydrograph  was  simulated  is 
also  shown. 

The  validity  of  the  computational  procedures  was  tested  by  reversing  the 
simulation  process.  The  previously  determined  inflows  (Figure  14)  were 
subjected  to  the  same  irrigation  demands  and  Irrigated  area  outflows  were 
simulated.  The  resulting  outflow  hydrograph  was  exactly  the  same  as  the 
historical  outflow  hydrograph  (Canyon  Ferry  inflow  hydrograph) .  Figure  15 
shows  a  portion  of  the  simulated  outflow  hydrograph.  It,  along  with  the 
inflow  hydrograph,  is  identical  to  the  hydrographs  of  Figure  14  (vertical 
scales  are  different) ,  thus  demonstrating  the  validity  of  the  computational 
procedures. 

Figure  15  provides  additional  insight  to  the  computational  procedure  for 
determining  irrigated  area  outflow.  The  high  peaked  hydrograph  (black  line) 
represents  the  "natural"  inflow  to  the  irrigated  area.  The  lowest  hydrograph 
(red  line)  represents  the  bypassed  flow  (the  "natural"  inflow  that  is  not 
diverted).  The  blue  lines  above  that  represent  groundwater  return  flows  from 
individual  months.  The  uppermost  blue  line  represents  the  irrigated  area 
outflow.  The  outflow  is  obtained  by  summing  the  bypassed  flow  and  the 
groundwater  return  flows  (see  Figure  16  for  clarification).  The  outflow 
(dashed  line)  in  September  is  obtained  by  summing  the  September  bypass  and  the 
return  flows  from  groundwater  recharges  during  September,  August,  and  previous 
months.  The  outflow  during  October  is  obtained  by  summing  the  October  bypass 
and  the  return  flows  from  the  groundwater  recharges  during  October,  September, 
August,  and  previous  months. 
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Irrigated   area   inflow    ("natural"    inflow)    as   determined 
from  known    irrigated  area  outflow. 


Figure  15.   Irrigated  area  outflow  as  determined  from  "natural" 
inflow.   Also  shown  are  the  groundwater  return  flows 
from  individual  months. 
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Figure  16.   Irrigated  area  outflow. 
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Figure  16.   Irrigated  area  outflow. 
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The  contribution  from  any  month's  recharge  diminishes  exponentially  with  time 
as  indicated  by  the  gradual  convergence  of  the  return  flow  lines  (Figure  15) . 

Although  the  aggregate  irrigated  area  inflows  and  outflows  of  Figure  16 
are  always  positive,  the  bypass  is  sometimes  negative.  The  negative  bypass  is 
a  computational  peculiarity  that  results  from  the  aggregation  of  the  irrigated 
areas.  Consider  Figure  17(a)  in  which  two  irrigated  areas,  A,  and  A_  are 
considered  separately.  The  symbols  I,  0,  D,  B,  and  R  represent  respectively 
inflow,  outflow,  diversion,  bypass  and  return  flow.  Repeated  application  of 


sv 


(b) 


Figure  17.  (a)  Nonaggregated  irrigated  areas  and  (b)  aggregated  irrigated 


Che  concinuicy  equation  in  the  downstream  direction  yields 


I-D^+B^ 
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°1  "  «1  ^  h 

0  =•  B2  +  R2 
It  is  noted  that  the  bypassed  flows  for  the  separate  areas  (B.  and  B.)  are  not 
negative.   It  is  not  possible  to  divert  more  water  from  a  point  of  diversion 
than  flows  to  that  point  of  diversion.   Combining  the  separate  continuity 
expressions  yields 

0  «  I  -  D^  +  R^  -  D^  +  Rj 

or  in  general  for  n  separate  irrigated  areas, 

n       n 
0«I-rD  +  rR 
1=1      i'l 
It  is  noted  here  that  the  total  diversions  can,  and  often  do,  exceed  the 

n 
Inflow  (Z  D.  >  I).   This  typically  occurs  in  the  late  summer  months,  and  it  is 
i=l  i 

possible  because  the  irrigated  area  inflow  is  augmented  within  the  irrigated 
area  by  return  flows  from  previous  months  irrigation  applications.  For 
example,  irrigated  area  A_  of  Figure  17(a)  can  divert  return  flows  from 
earlier  irrigation  applications  to  A.  as  well  as  the  current  month's  bypass  of 
A..  The  result  is  a  total  monthly  irrigated  area  diversion  (to  A  and  A„) 
that  exceeds  the  monthly  Irrigated  area  inflow. 

Now  consider  Figure  17(b),  a  schematic  representation  in  which  the  n 
Irrigated  areas  have  been  aggregated.  If  the  schematic's  undiverted  flow  is 
defined  as  the  bypass  B  (analogous  to  the  separate  irrigated  area  bypasses  of 
Figure  17(a)),  then  repeated  application  of  the  continuity  equation  in  the 
downstream  direction  yields 

1  =  Z  D,  +  B 

1=1  i 


0  =  B  +  Z  R. 

i  =  L  1 
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Combining  the  separate  continuity  expressions, 

n      n 

1=1  i   1=1  i 


The  equation  for  the  outflow  from  the  aggregate  area  is  identical  to  the 

equation  for  the  outflow  from  the  unaggregated  areas.   Since  this  study  is 

concerned  with  the  impact  of  the  total  irrigation  activity,  and  not  with 

impacts  within  the  aggregate  area,  the  aggregation  is  quite  acceptable. 

However,  as  was  pointed  in  the  previous  paragraph,  diversions  may  exceed 

n 
inflow  (I  D   >  I),  and  when  they  do,  the  so-called  aggregate  irrigated  area 

i  =  l  ^ 

bypass  of  Figure  17(b)  will  be  negative.  The  negative  bypass  is  thus  a 
peculiarity  that  results  from  the  aggregation  of  irrigated  areas.  It  does 
not,  however,  in  any  way  invalidate  the  analysis. 

With  reference  to  Figure  15,  it  can  be  seen  that  the  outflow  hydrograph 
is  less  variable  than  the  inflow  hydrograph.  Irrigation  season  outflows  have 
generally  been  reduced  by  irrigation  diversions  while  off-season  outflows  have 
been  increased  by  irrigation  return  flows.  Although  variable  from  year  to 
year,  the  hydrograph  "smoothing"  afforded  by  the  diversions  and  delayed  return 
is  always  quite  significant.  Early  summer  peaks  are  typically  reduced  by  25 
percent  or  so.  Late  summer  flows  may  be  reduced  by  70  percent  or  more. 
Winter  flows  may  be  increased  by  more  than  100  percent. 

The  hydrograph  smoothing  that  results  from  irrigation  activities  is  of 
major  importance.  Among  other  things,  it  greatly  reduces  the  possibility  of 
flooding  during  the  snow  melt  runoff  season.  It  also  benefits  power  producers 
in  that  streamflows  are  maintained  at  artificially  high  levels  during  the 
winter,  thus  permitting  increased  power  production  during  this  period  of  high 
demand.  If  the  hydrograph  smoothing  agent,  irrigation  activity,  is  changed, 
Che  extent  of  the  smoothing,  and  the  benefits  that  derive  therefrom  will,  of 
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course,  be  changed.   The  following  sections  investigate  the  effects  of  such 
changes. 

SIMULATIONS  OF  IRRIGATION  ACTIVITY  CHANGES 

The  irrigated  area  inflows  obtained  in  the  previous  section  were 
subjected  to  specified  irrigation  demands  and  the  irrigated  area  outflows  were 
simulated.  The  outflows  thus  obtained  were  then  input  to  the  Canyon  Ferry 
Reservoir  where  the  previously  discussed  reservoir  operation  rules  determined 
their  ultimate  disposal.  Two  general  types  of  irrigation  activity  change  were 
simulated:  changes  in  irrigation  efficiency,  and  changes  in  irrigated 
acreage. 

Simulations  of  Changes  in  Irrigation  Efficiency 

Irrigation  efficiencies,  which  have  been  increasing  with  expansions  in 
sprinkler  use,  reductions  in  conveyance  loss,  and  better  management,  are 
projected  to  increase  even  more  in  the  future  (1).  The  future  increases  will 
probably  be  gradual  and  perhaps  sporadic  as  economic  conditions  vary.  They 
will  almost  certainly  occur,  however,  as  irrigators  seek  to  improve  their 
profits. 

More  efficient  irrigation  operations  require  smaller  diversions  and 
produce  less  return  flow,  and  thereby  increase  irrigated  area  outflow 
variability.  Spring  and  early  summer  outflows  are  higher  because  of  the 
reduced  diversions  and  winter  outflows  are  lower  because  of  the  reduced 
groundwater  return  flows.   The  increased  variability  in  outflow  (Canyon  Ferry 
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Inflow)  impacts  reservoir  operations  as  they  relate  to  power  production,  flood 
control,  etc. 

The  effects  of  irrigation  efficiency  improvements  were  Investigated  bv 
simulating  the  Canyon  Ferry  system's  hydrologic  behavior  under  three  levels  of 
assumed  irrigation  efficiency.  All  other  system  variables  such  as  irrigated 
area  inflow,  crop  consumptive  needs,  etc.  were  held  constant,  so  that  the 
separate  effects  of  the  efficiency  changes  might  be  assessed.  Efficiencies  of 
19  (the  estimated  current  efficiency  (1)),  40,  and  60  percent  were  considered 
in  the  three  simulations.  It  is  recognized,  of  course,  that  the  higher 
efficiencies,  particularly  the  60  percent  efficiency,  might  never  be  attained. 
Simulations  of  these  efficiency  conditions  are  never-the-less  useful  in  that 
they  clearly  identify  the  character  and  relative  significance  of  the  efficien- 
cy changes. 

Simulation  results  were  summarized  in  a  series  of  graphs.  Portions  of 
these  graphs  (covering  the  period  from  1969  through  1973)  are  shown  in  the 
following  figures.  Each  of  the  graphs  includes  four  records:  the  historical 
record  (Bureau  of  Reclamation)  and  one  record  for  each  of  the  three  simulated 
irrigation  efficiencies. 

The  first  graph.  Figure  18,  shows  how  Canyon  Ferry  inflow  is  affected  by 
changes  in  irrigation  efficiency.  Because  the  simulated  inflow  hydrograph  for 
the  19  percent  irrigation  efficiency  (estimated  current  efficiency)  is 
coincident  with  the  historical  inflow  hydrograph,  only  3  hydrographs  are 
evident  in  the  figure.  The  reservoir  inflow  hydrograph  peaks  are  increased  as 
the  irrigation  efficiency  is  increased  and  correspondingly  less  water  is 
diverted.  As  will  be  shown  in  subsequent  figures,  the  increased  spring 
reservoir  inflows  necessitate  more  frequent  use  of  the  reservoir's  flood 
control  storage  as  well  as  reservoir  releases  that  more  frequently  exceed  the 
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below-dam  channel  capacity.  Figure  18  also  shows  that  winter  reservoir 
Inflows  decrease  with  increases  in  Irrigation  efficiency  and  the  resultant 
reductions  in  groundwater  return  flow.  Although  variable,  reductions  in 
Canyon  Ferry  winter  inflow  range  up  to  45  percent  and  more  when  irrigation 
efficiency  is  increased  from  19  to  60  percent.  As  will  be  shown  in  subsequent 
figures,  these  reduced  inflows  limit  energy  production  during  the  winter 
months. 

The  next  four  figures  (19,  20,  21,  and  22)  illustrate  how  Canyon  Ferry 
operations  are  impacted  by  the  changes  in  irrigation  efficiency  and  the 
resultant  changes  in  reservoir  in  flow.  In  these  figures,  the  simulated 
records  for  the  estimated  current  irrigation  efficiency  of  19  percent  (solid 
black  lines)  are  not  coincident  with  the  historical  records  (dashed  black 
lines)  .  These  differences  result  from  the  fact  that  the  operations  rules 
employed  in  the  simulations  deviate  somewhat  from  those  that  were  historically 
employed  in  the  operation  of  the  reservoir.  It  would  be  very  difficult, 
probably  impossible,  to  specify  a  set  of  reservoir  operations  rules  that 
would,  when  applied  in  a  simulation,  yield  reservoir  storages,  etc.  that  were 
identical  to  historical  records  of  storage,  etc.  The  historical  records 
reflect  a  rather  complex  and  evolving  set  of  flexible  reservoir  operations 
procedures  that  were  applied  in  the  daily  operation  of  the  reservoir.  The 
simulation  rules,  on  the  other  hand,  are  necessarily  simpler,  less  flexible, 
and  they  are  applied  on  a  monthly  basis.  These  simpler  rules  never-the-less 
produce  a  reasonable  simulation  of  historical  reservoir  operations  (compare 
solid  and  dashed  black  lines) .  The  reader  is  referred  to  the  section  titled 
"Reservoir  Operation  Procedures"  for  a  detailed  discussion  of  the  operations 
rules  and  the  rationale  behind  them. 
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Inflows  decrease  with  increases  in  irrigation  efficiency  and  the  resultant 
reductions  in  groundwater  return  flow.  Although  variable,  reductions  in 
Canyon  Ferry  winter  inflow  range  up  to  45  percent  and  more  when  irrigation 
efficiency  is  increased  from  19  to  60  percent.  As  will  be  shown  In  subsequent 
figures,  these  reduced  Inflows  limit  energy  production  during  the  winter 
months , 

The  next  four  figures  (19,  20,  21,  and  22)  illustrate  how  Canyon  Ferry 
operations  are  Impacted  by  the  changes  in  irrigation  efficiency  and  the 
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black  lines)  are  not  coincident  with  the  historical  records  (dashed  black 
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employed  in  the  simulations  deviate  somewhat  from  those  that  were  historically 
employed  in  the  operation  of  the  reservoir.  It  would  be  very  difficult, 
probably  impossible,  to  specify  a  set  of  reservoir  operations  rules  that 
would,  when  applied  in  a  simulation,  yield  reservoir  storages,  etc.  that  were 
Identical  to  historical  records  of  storage,  etc.  The  historical  records 
reflect  a  rather  complex  and  evolving  set  of  flexible  reservoir  operations 
procedures  that  were  applied  in  the  daily  operation  of  the  reservoir.  The 
simulation  rules,  on  the  other  hand,  are  necessarily  simpler,  less  flexible, 
and  they  are  applied  on  a  monthly  basis.  These  simpler  rules  never-the-less 
produce  a  reasonable  simulation  of  historical  reservoir  operations  (compare 
solid  and  dashed  black  lines)  .  The  reader  is  referred  to  the  section  titled 
"Reservoir  Operation  Procedures"  for  a  detailed  discussion  of  the  operations 
rules  and  the  rationale  behind  them. 
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Figure  19  shows  how  the  month-end  reservoir  storage  varies  with  changes 
in  irrigation  efficiency.  The  reservoir  operations  rules,  which  are  stated  on 
the  figure,  were  the  same  in  each  of  the  three  efficiency  simulations.  Thus 
the  plots  in  Figure  19  (and  in  Figures  20,  21  and  22)  illustrate  impacts  that 
are  solely  attributable  to  changes  in  irrigation  efficiency.  It  can  be  seen 
from  the  figure  that  month-end  storage  generally  becomes  more  variable  as 
irrigation  efficiency  increases  and  that  flood  control  storage  (above  upper 
horizontal  line)  must  be  more  frequently  and  more  extensively  used.  During 
the  period  from  1968  through  1984  ,  it  was  necessary  to  use  this  storage 
during  6,  11,  and  13  years  for  respective  irrigation  efficiency  simulations  of 
19,  40,  and  60  percent.  Figure  19  also  shows  that  reservoir  storage  falls 
below  that  level  at  which  head  limits  power  production  (middle  horizontal 
line)  more  frequently  as  irrigation  efficiency  increases. 

Figure  20  illustrates  how  the  discharge  below  the  dam  (total  reservoir 
release  to  the  Missouri  River)  varies  with  changes  in  irrigation  efficiency. 
It  can  be  seen  that  the  channel  capacity  is  exceeded  more  often  and  more 
extensively  as  irrigation  efficiency  increases.  The  reservoir  operations 
rules  permit  the  channel  capacity  to  be  exceeded  only  when  flood  control 
storage  (Figure  19)  is  exhausted.  During  the  period  from  1968  through  1984, 
it  was  necessary  to  exceed  channel  capacity  during  0,  8,  and  9  years  for 
respective  irrigation  efficiency  simulations  of  19,  40  and  60  percent. 

Figure  21  shows  how  power  turbine  discharge  is  affected  by  changes  in 
Irrigation  efficiency.  The  most  notable  effects  are  the  reductions  in  winter 
turbine  discharge  that  accompany  increases  in  irrigation  efficiency.  Although 
variable,  the  winter  discharge  may  be  reduced  as  much  as  50  percent  when  Che 
irrigation  efficiency  is  increased  from  19  to  60  percent.   This  is  primarily 
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Figure    19.      The  effects  of   changes   in    irrigation   efficiency  on  Canyon 
Ferry  month-end  storage. 


iflPACr  or  ACCRECATE  IRRlCATiON  r,C'i1U1    ON  CANYOX  FfRRr  RESERVOIR  OPERAIIONS 
THANNEL  DISCHARGE  BELOU  DAn 

OBSERVED  (BU-  RECl  RECORDS) 

CQUPUTEO  (BASED  ON  FOLLOWING  ASSUMED  RESERVOIR  OPERATION  RULES  •  RilLES  FIRST  APPLIED  IN  WATER  YEAR  1956  AFTER  RESERVOIR  INITIALLf  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTlNUOUSLr  SINCE  THEN  THOUGH  ONLV  A  PORTION  OF  THE  SIMULATION  MAV  BE  PLOTTED! 
I.  RESERVOIR  CAPACiT.r  (2050900  A-FT  AT  W  S  ELEV  3800  FT)  MUST  NOT  BE  EXCEEDED-  RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITY]  TO  ASSURE- 

^.  RESERVOIR  STORACE  SHOULD  NOT  EXCEED  1946624  A-FT  (W  S  ELEV  3797  FJ).    STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS.  RELEASE  WATER  IF  CHANNEL  CAPACITV  (15000  CFS) 
J-  RESERVOIR  STORACE  AT  MONTH-END  SHOULD  A^'PRQXIMATE  SPECiFlEO  MONTHLT  COAL-  RELEASE  WATER  TO  MEET  (APPROXIMATE)  COAl ■  RELEASES 

ARE  CONSTRAINED  Sf    A  MINIMUM  DOWNSTREAM  FLOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  i 5000  CFS- 
4-  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE)- 
CROP  CONSUMPTIVE  N££D=^CURRENT  CROP  CONSUMPTIVE  NEED 
IRRIGATION  EFFIClENCY-19/  X 

IRRIGATION  EFFL0LENCr>407  ■ 

lRRIGAT10>J  -;^FiClENC'.60:  X 

CHANNEL  CAPACITY  -  15000  CFS 
MINIMUM  FLOW  REQUIREMENT  -  2500  CFS 


The  effects  of  changes  in  irrigation  efficiency  on 
Canyon  Fsrry  releases  to  the  Missouri  River. 


IMPACT  OF  AGGREGATE  IRRIGATION  ACTIVITY  ON  CANYON  FERRY  RESERVOIR  OPERATIONS 
='OJER  TURBINE  DISCHARGE 

OBSERVED  (BU-  RECl   RECORDS) 

COriPUIEO  (BASED  ON  FOLLOWING  ASSjriED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  IN  HATER  YEAR  1956  AFTER  RESERVOIR  INITIALLY  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTINUOUSLY  SINCE  THEN  THOUGH  ONLY  A  PORTION  OF  THE  SIMULATION  MAY  BE  PLOTTED) 
,   RESERVOIR  CAPACITY  12050900  A-FT  AT  W  S  ELE/  3800  FTl  MUST  NOT  BE  EXCEEDED.  RELEASE  UATER  (RECARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITY]  TO  ASSURE 

<:■  RESERVOIR  STORAGE  SHOULD  NOT  EXCEED  1946621  A-FT  (W  S  ELEV  3'9r  FT).  STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS   RELEASE  WATER  IE  CHANNEL  CAPACITY  (15000  CFSl 
J.  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLY  GOAl -  RELEASE  WATER  TO  MEET  (APPROXIMATE)  COAL.  RELEASES 

ARE  CONSTRAINED  BY  A  MINIMUM  DOWNSTREAM  FLOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  15000  CFS. 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  [WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE). 
CROP  CONSUMPTIVE  NEEO=CURRENT  CROP  CONSUMPTIVE  NEED 
iRRlCATiON  EFFICIENCY-, 9«  M 

iRRlCAriON_EFFlClENCY-tOZ  Z 

IRRIGATION  £FFiCl£NC«  =  60-^  X 


Figure  21.   The  effects  of  changes  in  irrigation  efficiency  on 
Canyon  Ferry  power  turbine  discharge. 


because  of  the  reduced  winter  groundwater  returns  that  are  associated  with  the 
higher  irrigation  efficiency  (see  low  winter  reservoir  inflows  at  higher 
efficiencies  in  Figure  18).  The  increased  efficiency  has  little  effect  on 
turbine  releases  during  the  sunnner  months  in  most  years.  That  is  because  the 
reservoir  inflows  in  the  spring  are  usually  high  enough  (even  when  reduced  by 
the  increased  diversions  required  for  low  irrigation  efficiency  operations)  to 
satisfy  turbine  discharge  needs  during  these  months.  However,  during  very  dry 
years  such  as  1973  (see  low  spring  reservoir  inflows  in  Figure  18),  the  summer 
turbine  discharge  may  be  significantly  affected  by  changes  in  irrigation 
efficiency.  The  increased  spring  inflows  that  result  from  the  smaller  di- 
versions required  by  the  more  efficient  operation  permit  larger  summer  turbine 
releases  (Figure  21).  Such  increased  releases,  however,  draw  down  the  reser- 
voir and  this  action  generally  necessitates  a  reduction  in  turbine  releases 
during  the  following  winter. 

Energy  production  is  dependent  upon  turbine  operating  characteristics, 
turbine  discharge,  and  turbine  head  (and  thus  reservoir  storage).  These 
quantities  are  related  through  the  power  plant's  operating  characteristic 
curves  (Appendix  F) .  As  discussed  in  the  previous  paragraphs,  reservoir 
storage  and  turbine  discharge  are  affected  by  changes  in  irrigation  efficien- 
cy. It  follows  then  that  energy  production,  which  is  dependent  upon  reservoir 
storage  and  turbine  discharge,  is  also  affected  by  changes  in  irrigation 
efficiency.  Figure  22  Illustrates  how  energy  production  is  affected  by  such 
changes.  The  curves  of  this  Figure  are  quite  similar  to  those  of  Figure  21 
since  energy  production  is  strongly  dependent  upon  turbine  discharge.  It  can 
be  seen  from  Figure  22  that  winter  energy  production  generally  decreases  as 
irrigation  efficiency  increases.  This  is  primarily  a  consequence  of  the 
reduced  turbine  discharge  that  results  from  the  diminished  winter  groundwater 
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because  of  the  reduced  winter  groundwater  returns  that  are  associated  with  the 
higher  irrigation  efficiency  (see  low  winter  reservoir  inflows  at  higher 
efficiencies  in  Figure  18).  The  increased  efficiency  has  little  effect  on 
turbine  releases  during  the  summer  months  in  most  years.  That  is  because  the 
reservoir  inflows  in  the  spring  are  usually  high  enough  (even  when  reduced  by 
the  increased  diversions  required  for  low  irrigation  efficiency  operations)  to 
satisfy  turbine  discharge  needs  during  these  months.  However,  during  very  dry 
years  such  as  1973  (see  low  spring  reservoir  inflows  in  Figure  18),  the  summer 
turbine  discharge  may  be  significantly  affected  by  changes  in  irrigation 
efficiency.  The  increased  spring  inflows  that  result  from  the  smaller  di- 
versions required  by  the  more  efficient  operation  permit  larger  summer  turbine 
releases  (Figure  21).  Such  increased  releases,  however,  draw  down  the  reser- 
voir and  this  action  generally  necessitates  a  reduction  in  turbine  releases 
during  the  following  winter. 

Energy  production  is  dependent  upon  turbine  operating  characteristics, 
turbine  discharge,  and  turbine  head  (and  thus  reservoir  storage).  These 
quantities  are  related  through  the  power  plant's  operating  characteristic 
curves  (Appendix  F)  .  As  discussed  in  the  previous  paragraphs,  reservoir 
storage  and  turbine  discharge  are  affected  by  changes  in  irrigation  efficien- 
cy. It  follows  then  that  energy  production,  which  is  dependent  upon  reservoir 
storage  and  turbine  discharge,  is  also  affected  by  changes  in  irrigation 
efficiency.  Figure  22  Illustrates  how  energy  production  is  affected  by  such 
changes.  The  curves  of  this  Figure  are  quite  similar  to  those  of  Figure  21 
since  energy  production  is  strongly  dependent  upon  turbine  discharge.  It  can 
be  seen  from  Figure  22  that  winter  energy  production  generally  decreases  as 
irrigation  efficiency  increases.  This  is  primarily  a  consequence  of  the 
reduced  turbine  discharge  that  results  from  the  diminished  winter  groundwater 
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IMPACT  OF  ACCRECATf  IRRICATION  AC'Hirr  ON  CAitrON  TERR'  RESERVOIR  OPERATIONS 
ENERC  PRODUC'ION 

OBSERVED  !BU.  RECl   RECORDSI 

COriPUTED  (BASED  ON  FOLLOWING  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  IN  WATER  rEAR  1 5S6  AFTER  RESERVOIR  INlTiALLY  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTINUOUSLT  SINCE  THEN  THOUGH  ONL'  A  PORTION  OF  THE  SIMULATION  MAf  BE  PLOTTED) 
i.  RESERVOIR  CAPACITr  (^050900  A-FT  AT  W  S  ELEV  J800  FTl  MUST  NOT  BE  EXCEEDED   RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPAClTri  TO  ASSURE 

i-    RESERVOIR  STORAGE  SHOULD  HOT  EXCEED  1946624  A  FT  (W  S  ElE/  3'9?  FTl.  STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS   RElEASE  WATER  IF  CHANNEL  CAPAClTr  (15000  CSl  PERMITS. 
J.  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLT  GOAL.  RELEASE  WATER  TO  MEET  (APPROXIMATE!  GOAL   RELEASES 

ARE  CONSTRAINED  ST  A  MINIMUM  DOWNSTREAM  FlOW  REOUIREMFNT  OF  iSOO  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITf  OF  1 5000  CFS ■ 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THA'  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  A80VEI ■ 
CROP  CONSUMPTIVE  NEEO-rURRENT  CROP  CONSUMPTIVE  NEED 
IRRIGATION  EFFIClENCr-19;(  ^ 

IRRIGATION  FFFlClFNOr-iOZ  Z 


Figure  22.   The  effects  of  changes  in  irrigation  efficiency 
Canyon  Ferry  energy  production. 


return  flows  associated  with  the  increased  efficiencies.  It  can  also  be  seen 
that  during  a  very  dry  year  such  as  1973,  an  increased  level  of  summer  energy 
production  is  possible  at  the  higher  irrigation  efficiency  levels.  This  is 
primarily  due  to  the  increased  turbine  discharge  during  this  period,  and  as 
with  the  turbine  discharge,  it  generally  comes  at  the  expense  of  a  reduction 
in  energy  production  during  the  following  winter. 

Table  2  below  further  illustrates  how  energy  production  is  affected  by 
changes  in  irrigation  efficiency.  The  table  lists  the  average  monthly  energy 
production  in  millions  of  Kw-Hrs  for  the  period  from  1968  through  1984. 

TABLE  2 

Average  Energy  Production  (Kw-Hr  x  10  )  from  1968-84 

Month     Observed    Simulated  for     Simulated  For  Simulated  For 

Bureau  of   19%  Irrigation    40Z  Irrigation  60%  Irrigation 

Reclamation   Efficiency       Efficiency  Efficiency 

25.9  22.7 

28.0  23.7 

38.1  32.5 

34.6  31.7 

32.1  29.1 

32.1  30.0 

38.1  37.8 

42.5  42.8 

41.1  42.2 

44.8  44.8 

44.5  44.9 

38.1  38.1 

440.1  420.3 
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Oct 

35.4 

34.9 

Nov 

37.1 

39.3 

Dec 

41.7 

43.9 

Jan 

42,2 

43.4 

Feb 

36.8 

39.1 

Mar 

39.0 

40.6 

Apr 

37.3 

39.4 

May 

38.8 

41.6 

June 

39.0 

39.3 

July 

39.8 

40.9 

Aug 

35.3 

36.8 

Sept 

31.0 

34.0 

Annual 

453.3 

473.2 

return  flows  associated  with  the  increased  efficiencies.  It  can  also  be  seen 
that  during  a  very  dry  year  such  as  1973,  an  increased  level  of  summer  energy 
production  is  possible  at  the  higher  irrigation  efficiency  levels.  This  is 
primarily  due  to  the  increased  turbine  discharge  during  this  period,  and  as 
with  the  turbine  discharge,  it  generally  comes  at  the  expense  of  a  reduction 
in  energy  production  during  the  following  winter. 

Table  2  below  further  illustrates  how  energy  production  is  affected  by 
changes  in  irrigation  efficiency.  The  table  lists  the  average  monthly  energy 
production  in  millions  of  Kw-Hrs  for  the  period  from  1968  through  1984. 

TABLE  2 

Average  Energy  Production  (Kw-Hr  x  10  )  from  1968-84 

Month     Observed    Simulated  for     Simulated  For  Simulated  For 

Bureau  of   19%  Irrigation    40%  Irrigation  603:  Irrigation 

Reclamation   Efficiency      Efficiency  Efficiency 

25.9  22.7 

28.0  23.7 

38.1  32.5 

34.6  31.7 

32.1  29.1 

32.1  30.0 

38.1  37.8 

42.5  42.8 

41.1  42.2 

44.8  44.8 

44.5  44.9 

38.1  38.1 

440.1  420.3 
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Oct 

35.4 

34.9 

Nov 

37.1 

39.3 

Dec 

41.7 

43.9 

Jan 

42.2 

43.4 

Feb 

36.8 

39.1 

Mar 

39.0 

40.6 

Apr 

37.3 

39.4 

May 

38.8 

41.6 

June 

39.0 

39.3 

July 

39.8 

40.9 

Aug 

35.3 

36.8 

Sept 

31.0 

34.0 

Annual 

453.3 

473.2 

Ic  can  be  seen  chat  che  average  annual  energy  production  (bottom  of  table)  is 
reduced  as  irrigation  efficiency  is  increased.  More  importantly,  the  time 
distribution  of  monthly  energy  production  is  significantly  changed  by  the 
irrigation  efficiency  changes.  In  general,  winter  energy  production  is 
reduced  and  summer  energy  production  is  increased  as  irrigation  efficiency  is 
increased. 

As  evidenced  by  Figures  18  through  22  and  by  Table  2,  it  must  be  conclud- 
ed that  increases  in  irrigation  efficiency  can  have  a  rather  significant 
impact  on  Canyon  Ferry  inflow  and  thereby  Canyon  Ferry  operations. 

Simulations  of  Changes  in  Irrigated  Acreage 

There  is  a  demand  for  more  extensive  use  of  the  surface  water  resource  in 
the  Canyon  Ferry  drainage  basin.  Applications  for  such  use  (primarily  for 
irrigation)  are  periodically  received  and  reviewed  by  the  state's  Water  Rights 
Bureau. 

An  increase  in  irrigated  acreage  requires  an  Increase  in  the  irrigation 
diversion  (assuming  no  change  in  irrigation  efficiency).  A  portion  of  the 
diversion  increase  is  consumed  within  the  newly  irrigated  area  and  a  portion 
of  it  returns  to  the  stream  system  via  surface  runoff  and  groundwater  flow. 
The  consumed  portion  is  of  course  lost  to  the  system  and  total  irrigated  area 
outflow  is  correspondingly  reduced.  The  surface  runoff  component  of  the 
diversion  increase  returns  quickly  to  the  stream  resulting  in  only  a  minimal 
change  in  irrigated  area  outflow  timing  and  no  change  in  outflow  volume.  The 
groundwater  return  component  of  the  diversion  increase  returns  very  slowly  to 
the  stream  resulting  in  a  more  significant  change  in  outflow  timing  but  no 
change  in  outflow  volume.   Outflows  are  reduced  during  the  irrigation  season 


71 


by  the  diversion  of  this  water  and  increased  during  the  nonirrigation  season 
as  it  gradually  returns  to  the  stream  system.  Reservoir  operations  (energy 
production,  flood  control,  etc.)  are  in  turn  impacted  by  these  changes  in  the 
volume  and  timing  of  irrigated  area  outflow  (Canyon  Ferry  inflow). 

The  effects  of  increases  in  the  irrigated  area  were  investigated  by 
simulating  the  Canyon  Ferry  system's  hydrologic  behavior  for  three  assumed 
irrigated  acreage  totals.  All  other  system  variables  such  as  irrigated  area 
inflow,  irrigation  efficiency,  etc.  were  held  constant  so  that  the  separate 
effects  of  the  acreage  changes  (and  proportionate  consumptive  use  changes) 
might  be  assessed.  Irrigated  acreages  of  658,900  (estimated  current  acreage 
(1)),  691,845,  and  724,790  acres  were  considered  in  the  simulations.  The 
latter  two  acreages  are  respectively  105  and  110  percent  of  the  first  (cur- 
rent) acreage. 

As  before,  simulation  results  were  summarized  in  a  series  of  graphs. 
Portions  of  these  graphs  (covering  the  period  from  1969  through  1973)  are 
shown  in  the  following  figures.  Each  of  the  graphs  include  four  records:  the 
historical  record  (Bureau  of  Reclamation)  and  one  record  for  each  of  the  three 
simulated  acreage  totals. 

The  first  graph.  Figure  23,  shows  how  the  Canyon  Ferry  inflow  is  affected 
by  changes  in  irrigated  area.  Because  the  simulated  hydrograph  for  the  es- 
timated current  acreage  is  coincident  with  the  historical  hydrograph,  only 
three  hydrographs  are  evident  in  the  figure.  Although  not  terribly  obvious 
from  the  figure,  the  reservoir  inflow  volumes  (areas  under  the  various 
hydrographs)  are  decreased  as  the  irrigated  acreage  is  increased  and 
consumptive  use  is  increased.  Somewhat  more  obvious  are  the  effects  that  the 
acreage  changes  have  on  the  timing  of  reservoir  inflows.  Irrigation  season 
inflows  are  reduced  by  the  increased  diversions  required  for  the  increased 
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iMPAcr  OF  accrecate:  irrigation  Acrwirr  on  canvon  FERRr  resfrvoir  operations 

rANrON  FERRY  INFlOJ 

OBSFR/FD  (BU.  RECL.  RECORDS) 

COMPUTED  (BASED  ON  FOL^OHING  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  iN  WATER  YEAR  1956  AFTER  RESERVOIR  INITIALLY  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTINUOUSLY  SINCE  THEN  THOUGH  ONLY  A  PORTION  OF  'HE  SiMULATION  MAY  BE  PLOrrEDI 
1-  RESERVOIR  CAPACITY  (2050900  A-FT  AT  U  S  ElEV  3800  FTI  MUST  NOT  BE  EXCEEDED.  RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITY)  TO  ASSURE. 

i-    RESERVOIR  SrORACE  SHOULD  NOT  EXCEED  1946624  A-FT  (W  S  ELEV  3?97  FT).  STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS.  RELEASE  WATER  IF  CHANNEL  CAPACITY  (15000  CFS)  PERMITS- 
J-  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLY  GOAl •  RELEASE  WATER  TO  MEET  (APPROXIMATE)  COA^ .  RELEASES 

ARE  CONSTRAINED  BY  A  11N1MUM  DOWNSTREAM  FLOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  15000  CFS 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE) • 
IRRIGATION  EFFICIENCY,  i  9/ 

iRR-  DIVERSlON-rORRENT  DIVERSION  X 

iRR.  DLVERSION'IOSZ  OF  CURRENT  DIVERSION    Z 

iRR-  DIVERS101-1I0-!  :.•■  ■■  .!-  =  ■:. 1"  ■jv<t's,z-t        y 


Figure    23.      The   effects   of   changes    in 
inflow. 


igate(3  area  on  Canyon  Ferry 


acreage  and  off-season  inflows  are  Increased  by  the  groundwater  return  flows 
from  those  increased  diversions. 

The  next  four  figure  (24,  25,  26,  and  27)  illustrate  how  Canyon  Ferry 
operations  are  impacted  by  the  changes  in  irrigated  acreage  and  the  resultant 
changes  in  reservoir  inflow.  In  these  figures,  the  simulated  records  for  the 
estimated  current  irrigated  acreage  (solid  black  lines)  are  not  coincident 
with  the  historical  records  (dashed  black  lines).  As  discussed  in  detail  in 
the  section  on  the  irrigation  efficiency  simulations,  these  differences  result 
from  the  fact  that  the  operations  rules  employed  in  the  simulations  necessari- 
ly deviate  somewhat  from  those  that  were  historically  employed  in  the  opera- 
tion of  the  reservoir.  The  simulation  rules  never-the-less  produce  a 
reasonable  simulation  of  historical  reservoir  operations  (compare  solid  and 
dashed  black  lines).  More  importantly,  the  effects  of  changes  in  irrigated 
area  are  more  apparent  when  viewed  in  the  framework  of  a  fixed  set  of 
reservoir  operations  rules  (compare  red,  blue,  and  solid  black  lines). 

Figure  24  shows  how  the  month-end  reservoir  storage  is  affected  by 
changes  in  irrigated  acreage.  Reservoir  storage  generally  decreases  in  the 
summer  months  and  increases  in  the  winter  months  in  response  to  the  changes  in 
reservoir  inflows  that  result  from  changes  in  irrigated  acreage.  The 
decreased  summer  storages  associated  with  the  increased  acreages  require  less 
frequent  and  less  extensive  use  of  flood  control  storage  (above  upper 
horizontal  line).  During  the  period  from  1968  through  1984,  it  was  necessary 
to  use  this  storage  during  6,  4,  and  1  years  in  respective  simulations  of  100, 
105,  and  110  percent  of  the  current  estimated  irrigated  acreage.  Figure  24 
also  shows  that  reservoir  storage  falls  below  that  level  at  which  head  limits 
power  production  (middle  horizontal  line)  less  frequently  as  the  irrigated 
area  increases. 
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acreage  and  off-season  inflows  are  increased  by  the  groundwater  return  flows 
from  those  Increased  diversions. 

The  next  four  figure  (24,  25,  26,  and  27)  illustrate  how  Canyon  Ferry 
operations  are  impacted  by  the  changes  in  irrigated  acreage  and  the  resultant 
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from  the  fact  that  the  operations  rules  employed  in  the  simulations  necessari- 
ly deviate  somewhat  from  those  that  were  historically  employed  in  the  opera- 
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also  shows  that  reservoir  storage  falls  below  that  level  at  which  head  limits 
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IMPACT  OF  ACCRECATE  IRRICAIION  ACTIVUI'  ON  CANYON  FERRY  RESERVOIR  OPERATIONS 
•10NTH-END  RESERVOIR  S'ORAGE 

OBSERVED  (SU-  RECL   RECORDS)         

COnPUTED  (BASED  OH  FOLLOUING  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  ^IRST  APPLIED  IN  WATER  TEAR  1956  AFTER  RESERVOIR  iNITlALLT  FIlLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTINUOUSLT  SINCE  THEN  THOUGH  ONLY  A  PORTION  OF  THE  SIMULATION  MAY  BE  PLOTTEOI 
;.  RESERVOIR  CAPACITY  (2050900  A-FT  AT  J  S  ELEV  3800  FT)  MUST  NOT  BE  EXCEEDED-  RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITTI  TO  ASSURE. 

,;.  RESERVOIR  STORAGE  SHOULD  NOT  EXCEED  1946624  A-FT  (W  S  ELEV  3797  FT).  STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS.  RELEASE  WATER  IF  CHANNEL  CAPACITV  (15000  CFS)  PERMITS. 
J   RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLT  COAL.  RELEASE  WATER  TO  MEET  'APPROXIMATE)  GOAL   RELEASES 

ARE  CONSTRAINED  BY  A  MINIMUM  DOWNSTREAM  FLOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  15000  CFS. 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE). 
IRRIGATION  EFFICIENCY'  i 9Z 
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MAXIMUM  RESERVOIR  STORAGE  (EXCLUDING  FLOOD  CONTROL  STORAGE)  »  1945624  A-FT. 
HEAD  LIMITS  POWER  GENERATION  BELOW  STORAGE  1497289  A-FT. 
"OWERPLANT  INOPERABLE  BELOW  STORAGE  434537  A-FT. 


AiV     fitt-  J«l' 


Figure  24.   The  effects  of  changes  in  irrigated  area  on  Canyon  Ferry 
month-end  storage. 


Figure  25  illustrates  how  the  discharge  below  the  dam  (total  reservoir 
release  to  the  Missouri  River)  varies  with  changes  in  irrigated  area.  Al- 
though never  exceeded  during  the  1969  through  1973  period,  the  channel  capaci- 
ty is  more  nearly  approached  when  the  smaller  areas  are  irrigated.  This  is 
because  a  smaller  portion  of  the  spring  irrigated  area  inflow  is  diverted  to 
the  smaller  irrigated  areas.  Figure  25  also  reveals  that  during  the  late 
summer  months,  the  reservoir  releases  to  the  river  are  smaller  when  the  larger 
areas  are  irrigated,  and  occasionally  may  fall  to  the  minimum  channel  flow 
requirement  of  2500  cfs  (lower  horizontal  line) .  This  is  because  the  in- 
creased diversions  required  for  the  larger  areas  diminish  reservoir  inflow 
significantly  during  the  low  flow  late  summer  months  causing  reservoir  re- 
leases to  be  reduced  via  the  reservoir  operations  rules.  During  the  fall  and 
winter,  releases  to  the  river  are  increased  marginally  when  the  larger  areas 
are  irrigated. 

Figure  26  shows  how  the  power  turbine  discharge  is  affected  by  changes  in 
Irrigated  areas.  It  is  little  affected  during  late  spring  and  early  summer. 
Snowmelt  runoff  peaks  at  this  time,  and  turbine  discharge  is  limited  by 
turbine  capacity  rather  than  reservoir  inflow  and  storage  regardless  of  the 
acreage  irrigated.  During  late  summer  and  early  fall,  the  turbine  discharge 
generally  decreases  with  an  increase  in  irrigated  areas  (and  resultant 
decrease  in  reservoir  inflow).  During  the  late  fall,  winter,  and  early 
spring,  there  is  a  net  increase  in  turbine  discharge  with  an  increase  in  the 
irrigated  area.  This  increased  discharge  is  made  possible  by  the  increased 
groundwater  return  flows  from  the  larger  areas.  The  off-season  increases  in 
turbine  discharge  occur  in  the  late  fall  of  each  year  and  in  the  early  spring 
of  some  years  (1970  and  1973  of  Figure  26).  The  time  distribution  of  these 
increased  releases  are  a  consequence  of  the  reservoir  operations  rules, 
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iNPACr  OF  ACCRECAIE  IRRlCAflON  ACTIKirv  ON  CANrON  FERRY  RESERVOIR  OPERATIONS 
CHANNEL  DISCHARGE  BELOW  DAM 

OBSERVED  (BU-  RECL  RECORDS) 

COnPUrED  (BASED  ON  FOLLOWING  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  IN  WATER  VFAR  1956  AFTER  RESERVOIR  INITIALLY  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONTINUOUSLY  SINCE  THEN  THOUGH  ONLY  A  PORTION  OF  THE  SIMULATION  MAY  BE  PLOTTED) 
..  RESERVOIR  CAPACITY  (2050900  A-FT  AT  W  S  ElEV  3800  FT)     MUST  NOT  BE  EXCEEDED-  RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITY)  TO  ASSURE- 

2-  RESERVOIR  STORAGE  SHOULD  NOT  EXCEED  1946624  A-FT  (W  S  ELEV  3797  FT).  STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS-  RELEASE  WATER  IF  CHANNEL  CAPACITY  (15000  CFS)  PERMITS- 
J-  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLY  GOAL   RELEASE  WATER  TO  MEET  (APPROXIMATE)  GOAL-  RELEASES 

ARE  CONSTRAINED  BY  A  MINIMUM  DOWNSTREAM  FuOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  15000  CFS- 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE) - 
IRRIGATION  EFFICIENCY^  ,9/ 

iRR.  DIVERSION=CURRENT  DIVERSION  X 

iRR.  0I.VERS10N.i05Z  OF  CURRENT  DIVERSION     Z 
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The   effects  of   changes    in   irrigated   area 
releases   to   the  Missouri   River. 


Canyon  Ferry 


IMPACT  OF  ACGRECATE  IRRlGAriON  AC'I'/irr  ON  CANfON  FERRY  RESERVOIR  OPERAriONS 
"OUER  TURBINE  DiSCHARCE 

OBSERVED  (BU-  RECL •  RECORDS)        

COMPUTED  (BASED  ON  FOLLOUiNG  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  IN  WATER  TEAR  1956  AFTER  RESERVOIR  INITiALLT  FKLED  -  RESERVOIR 
OPERATION  SIMULATED  CONT 1 NUOUSLT  SINCE  THEN  THOUGH  ONLT  A  PORTION  OF  THE  SIMULATION  MAT  BE  PLOTTED) 
..  RESERVOIR  CAPACITr  (2050900  A-FT  AT  J  S  ELEV  3800  FT)  MOST  NOT  BE  EXCEEDED   RELEASE  WATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITf)  TO  ASSURE. 

Z      RESERVOIR  STORAGE  SHOULD  NOT  EXCEED  1946624  A-FT  (W  S  ELE^  3?9'  FT)   STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS   RELEASE  WATER  IF  CHANNEL  CAPACITT  (ISOOO  C^S)  PERMITS- 
J.  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHL)"  COAL-  RELEASE  WATER  TO  MEET  (APPROXIMATE)  COAL-  RELEASES 

ARE  CONSTRAINED  BT  A  MINIMUM  DOWNSTREAM  FLOW  REOUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  1 SOOO  CFS • 
■>-  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE!  - 
IRRIGATION  EFFICIENCY"  i 9Z 
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Figure  26.   The  effects  of  changes  in  irrigated  area  on  Canyon  Ferry 
power  turbine  discharge. 


particularly  the  requirement  that  specified  month-end  storage  goals  be 
approximated.  The  increased  turbine  discharge  could  be  released  at  any  time 
during  the  winter  with  an  appropriate  change  in  the  operations  rules.  Release 
during  the  high  energy  demand  winter  months  would  of  course  be  advantageous. 

Figure  27  demonstrates  how  energy  production  is  affected  by  changes  in 
irrigated  area.  Energy  production  is  strongly,  though  not  totally,  dependent 
upon  turbine  discharge,  and,  as  a  consequence,  the  factors  that  affect  turbine 
discharge  similarly  affect  energy  production.  Energy  production  is  largely 
unaffected  by  changes  in  irrigated  area  during  the  late  spring  and  early  fall 
when  snowmelt  runoff  occurs  and  reservoir  inflow  is  high.  It  generally 
decreases  in  the  late  summer  and  early  fall  when  the  increased  diversions 
required  for  the  larger  areas  diminish  reservoir  inflows,  and  Increases  during 
the  late  fall,  winter  and  early  spring  when  reservoir  inflows  are  increased  by 
the  groundwater  return  flows  derived  from  the  earlier  increased  diversions. 

Table  3  further  illustrates  the  effects  that  changes  in  irrigated  acreage 
have  on  energy  production.  The  Table  lists  average  monthly  energy  and  produc- 
tion in  millions  of  Kw-Hrs.  It  is  based  on  simulations  for  the  period  from 
1968  through  1984. 
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irtPACl  OF  ACGRECATE  IRRlCAnON  ACWUr  ON  CANTON  FERRY  RESERVOIR  OPERATIONS 
FNEROr  PRODUCTION 

OBSERVED  (BU  RECL   RECORDS!        

COMPUTED  (BASED  ON  FOLLOUINC  ASSUMED  RESERVOIR  OPERATION  RULES  -  RULES  FIRST  APPLIED  IN  WATER  YEAR  1956  AFTER  RESERVOIR  INITiALLf  FILLED  -  RESERVOIR 
OPERATION  SIMULATED  CONT 1 NUOUSLr  SINCE  THEN  IHOUCH  ONLY  A  PORTION  OF  THE  SIMULATION  MAY  BE  PLOTTED) 
..  RESERVOIR  CAPACiTY  (2050900  A-FT  AT  J  S  ELEV  J800  FTl  MUST  NOT  BE  EXCEEDED-  RELEASE  UATER  (REGARDLESS  OF  DOWNSTREAM  CHANNEL  CAPACITY!  TO  ASSURE. 

I      RESERVOIR  STORAGE  SHOULD  NOT  EXCEED  1946624  A-FT  (U  S  ELEV  3791'  FTl   STORAGE  RESERVED  FOR  FLOOD  CONTROL  ABOVE  THIS   RELEASE  WATER  IF  CHANNEL  CAPACITY  (15000  CFSI  PERMITS- 
J-  RESERVOIR  STORAGE  AT  MONTH-END  SHOULD  APPROXIMATE  SPECIFIED  MONTHLY  COAL   RELEASE  WATER  TO  MEET  (APPROXIMATE)  COAL-  RElEASES 

ARE  CONSTRAINED  BY  A  MINIMUM  DOWNSTREAM  FLOW  REQUIREMENT  OF  2500  CFS  AND  A  DOWNSTREAM  CHANNEL  CAPACITY  OF  15000  CFS ■ 
4.  RESERVOIR  RELEASES  SHOULD  BE  MADE  SUCH  THAT  POWER  TURBINE  DISCHARGE  IS  MAXIMIZED  (WITHIN  THE  RESERVOIR  STORAGE  AND  CHANNEL  FLOW  CONSTRAINTS  SPECIFIED  ABOVE) - 
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Figure  27. 


The  effects  of  changes  in  irrigated  area  on  Canyon  Ferry 
energy  production. 


TABLE  3 

Average  Energy  Production  (Kw-Hr  x  10  )  from  1968-1985 

Month     Observed       Simulated  for       Simulated  for     Simulated  for 

Bureau         Current  Irr.        105^  of  Current    110%  of  Current 
Reclamation    Acreage Irr.  Acreage Irr.  Acreage 

35.4  36.1 

39.5  40.1 

44.0  43.1 
43.5  43.9 

39.3  39.9 
41.7  43.2 
39.7  40.9 
41.5  41.5 

39.1  38.6 
40.1  39.1 

34.4  29.9 

33.5  32.1 
471.8             468.5 

It  can  be  seen  that  the  average  annual  energy  production  (bottom  of  table)  is 
decreased  only  slightly  when  the  acreage  Irrigated  is  increased  (from  473.2  x 
10  Kw-Hrs  for  the  current  acreage  to  468.5  x  10  Kw-Hrs  for  110  percent  of 
the  current  acreage) .  It  can  also  be  seen  that  the  time  distribution  of 
energy  production  is  changed  somewhat  when  the  irrigated  acreage  is  changed. 
A  comparison  of  the  simulation  results  reveals  that  during  the  October  through 
April  off  season,  an  average  of  280.6  x  10  Kw-Hrs  was  produced  when  the 
current  acreage  was  irrigated  and  an  average  of  287.2  x  10  Kw-Hrs  was 
produced  when  110  percent  of  the  current  acreage  was  irrigated. 
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Oct 

35.4 

34.9 

Nov 

37.1 

39.3 

Dec 

41.7 

43.9 

Jan 

42.2 

43.4 

Feb 

36.8 

39.1 

Har 

39.0 

40.6 

Apr 

37.3 

39.4 

May 

38.8 

41.6 

June 

39.0 

39.3 

July 

39.8 

40.9 

Aug 

35.3 

36.8 

Sept 

31.0 

34.0 

Annual 

453.2 

473.2 
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33.5  32.1 
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A  comparison  of  the  simulation  results  reveals  that  during  the  October  through 
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Oct 

35.4 

34.9 

Nov 

37.1 

39.3 

Dec 

41.7 

43.9 

Jan 

42.2 

43.4 

Feb 

36.8 

39.1 

Mar 

39.0 

40.6 

Apr 

37.3 

39.4 

May 

38.8 

41.6 

June 

39.0 

39.3 

July 

39.8 

40.9 

Aug 

35.3 

36.8 

SepC 

31.0 

34.0 

Annual 

453.2 

473.2 

As  evidenced  by  Figures  23  through  27  and  by  Table  3,  it  is  concluded 
that  increases  in  irrigated  acreage  have  a  minimal  and  not  necessarily  nega- 
tive impact  on  Canyon  Ferry  operations. 


82 


SUMMARY 
A  methodology  was  developed  to  simulate  the  hydrologlc  behavior  of  the 
Canyon  Ferry  system.   The  methodology  was  utilized  to  investigate  the  impacts 
of  changes  in  irrigation  efficiency  and  of  changes  in  irrigated  acreage. 

Improvements  in  irrigation  efficiency,  though  advantageous  to  agricul- 
tural interests,  were  shown  to  be  generally  disadvantageous  with  regard  to 
reservoir  operations.  Most  importantly  the  potential  for  spring  flooding  is 
increased  and  the  potential  for  energy  production,  particularly  winter  produc- 
tion, is  decreased  as  irrigation  efficiency  is  increased. 

Increases  in  irrigated  acreage,  again  advantageous  to  agricultural 
interests,  were  shown  to  have  minimal  and  mixed  impacts  with  regard  to  reser- 
voir operations.  On  the  negative  side,  annual  energy  production  potential  is 
reduced  slightly  as  irrigated  acreage  is  increased.  On  the  positive  side, 
spring  flood  potential  is  reduced  slightly  and  winter  energy  production 
potential  is  increased  slightly  as  the  irrigated  area  is  increased. 

As  detailed  in  earlier  sections,  numerous  assumptions  and  idealizations 
were  incorporated  in  the  methodology.  Though  theoretically  unappealing,  they 
are  a  practical  necessity.  It  would  be  impossible  to  simulate  the  behavior  of 
a  system  as  large  and  as  diverse  as  the  Canyon  Ferry  system  without  making 
numerous  assumptions.  Although  the  assumptions  employed  in  this  study  are 
believed  to  be  appropriate,  it  is  recognized  that  they  are  somewhat  subjective 
and  that  simulation  results  are  dependent  upon  them.  Never-the-less,  it  is 
believed  that  the  impacts  of  irrigation  activity  changes  identified  in  the 
study,  though  perhaps  not  precise,  are  never-the-less  quantitatively  reason- 
able and  qualitatively  valid.  That  is  to  say,  the  nature  of  the  impacts  would 
not  be  changed  by  reasonable  changes  in  the  assumptions. 
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Most  of  the  effort  expended  to  date  in  this  investigation  has  been 
directed  toward  the  development  of  a  methodology  for  assessing  the  impacts  of 
changes  in  irrigation  activity.  Many  other  things  could  and  should  be  done  to 
extend  the  usefulness  of  the  study.  The  effects  of  combined  changes  in 
irrigation  efficiency  and  irrigated  acreage  should  be  investigated  since  such 
combined  changes  are  likely  to  occur.  Some  combination  or  combinations  of 
these  two  changes  may  be  most  beneficial  to  the  two  major  water  users: 
agriculture  and  hydropower.  The  sensitivity  of  the  simulation  results  to 
changes  in  assumptions  relating  to  aquifer  properties,  irrigation  operations, 
reservoir  operations,  etc.  should  be  examined  as  well.  The  impacts  of  changes 
in  irrigation  activities  on  downstream  energy  producers  should  also  be  con- 
sidered since  these  downstream  users  are  very  much  dependent  upon  Canyon  Ferry 
releases  for  their  operation.  Optimization  studies  could  also  be  conducted  to 
identify  how  available  water  might  be  best  used  to  maximize  the  combined 
economic  benefits  to  the  agricultural  and  hydropower  users.  The  methodology, 
which  is  quite  general  might  also  be  applied  at  other  western  sites  where 
agriculture  and  hydropower  compete  for  a  limited  water  resource. 
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APPENDIX  A 
Table  A-1  lists  irrigation  figures  for  the  Canyon  Ferry  system.   All 
values  are  from  the  Water  Conservation  and  Salvage  Report  for  Montana  (1). 

TABLE  A-1 

Upper  (above  Canyon  Ferry)  Missouri  River 

Irrigation  Figures 


County 

Irrigated 

Crop  Irr. 

Irrigation 

Water 

Water 

Area 

Consumption 

Diversion 

(lO^Acres) 

(10^  A-ft) 

(10^  A-ft) 

Beaverhead 

329.0 

323.6 

1712.0 

Broadwater 

49.7 

60.3 

305.0 

Deer  Lodge 

8.6 

7.7 

45. 5 

Gallatin 

93.0 

78.3 

373.1 

Jefferson 

27.0 

27.5 

131.1 

Madison 

1A0.5 

144.3 

881.0 

Meagher 

6.5 

6.2 

42.0 

Silverbow 

4.6 

3.4 

13.7 

Totals  658.9  651.3  3503.4 

The  fraction  of  the  Meagher  County  Missouri  River  Basin  irrigated  acreage  that 

drains  Into  the  Canyon  Ferry  Reservoir  was  estimated  from  the  Atlas  of  Water 

Resources  in  Montana  (4)  to  be  1/7.   On  the  basis  of  the  figures  above  the 

irrigation  efficiency  within  the  basin  was  estimated  to  be 

Irrigation  Efficiency  -        x  lOOZ  s  19Z 
3503.9 


86 


Records  of  irrigation  diversions  at  several  sites  in  the  Upper  Missouri 
River  Basin  were  examined  to  determine  diversion  timing.  The  diversion 
records  are  summarized  in  Table  A-2  below. 

TABLE  A-2 

Average  Monthly  Irrigation  Diversions 
May  Through  September 


Helena  Valley 
From  Canyon  Ferry 
(1959-81) 

East  Bench 
Full  Service 
Near  Dillon 
(1966-73) 

East  Bench 
Supplemental 
Near  Dillon 
(1965-73) 

Volume 
Diverted 

Percent 

Volume 
Diverted 

Percent 

Volume 
Diverted 

Percent 

Month 

(acre- 

ft) 

of 

Total 

.   (acre-ft) 

of  total 

(acre-ft 

of  total 

May 

11620 

17.7 

11293 

13.4 

23085 

11.7 

June 

149A0 

22.8 

18933 

22.4 

55566 

28.1 

July 

17064 

26.0 

22017 

26.1 

45072 

22.8 

Aug 

15288 

23.3 

20104 

23.8 

47493 

24.0 

Sept 

6658 

10.2 

12077 

14.3 

26486 

13.4 

Totals: 

65570 

100.0 

84424 

100.0 

197702 

100.0 

These  records  were  assumed  to  be  representative  of  irrigation  timing  through- 
out the  basin  (records  are  generally  not  available  for  most  diversions).  On 
the  basis  of  these  records,  it  was  assumed  that  25  percent  of  the  total  annual 
diversion  occurs  in  each  of  the  three  months  of  June,  July,  and  August,  while 
10  percent  occurs  in  each  of  the  months  of  May  and  September,  and  2  1/2 
percent  occurs  in  April  and  October  (although  limited  diversions  are  sometimes 
made  in  April  and  October,  incomplete  records  made  it  necessary  to  estimate 
the  diversions  during  these  months) . 
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Table  A-3  lists  the  total  basln-wlde  diversions  by  month.  It  also  lists 
the  crop  consumptive  use  by  month  (assuming  an  irrigation  efficiency  of  19 
percent) . 

TABLE  A-3 
Monthly  Irrigation  Figures  for  Canyon  Ferry  Catchment 


Month 


Percent  of 
Annual  Diversion 


Diversion 
(Acre-ft) 


Crop  Irr. 

Water 
Consumption 
(Acre-ft) 


Apr 

2.5 

May 

10.0 

June 

25.0 

July 

25.0 

Aug 

25.0 

SepC 

10.0 

Oct 

2.5 

Annual 

Totals 

100 

87,585 
350, 3A0 
875,850 
875,850 
875,850 
350,340 

87,585 

3,503,400 


16,282 

65,130 

162,825 

162,825 

162,825 

65,130 

16,282 

651,300 
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APPENDIX  B 

Aquifer  transmissivity  values  determined  by  USGS  well  tests  in  the 
Gallatin  Valley  (5)  are  shown  in  Tables  B-1  and  B-2.  These  values  were  used 
to  estimate  the  average  tansmissivity  in  the  Canyon  Ferry  drainage. 

TABLE  B-1 

Transmissivity  Values  for 
Wells  Located  on  the  Alluvial  Plain 


Township 


Rang( 


Section 


Location  in 
Sections 


Transmissivity 
(gal/day/ft) 


2S 
3S 

IN 
IS 
IS 
IS 
IS 
IS 
IS 
IS 
IS 
2S 
25 
2S 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 


AE 
4E 
4E 
4E 
4E 
4E 
4E 
4E 
5E 
5E 
5E 
4E 
4E 
4E 
4E 
4E 
4E 
4E 
4E 
4E 
3E 
3E 
3E 
3E 


26 

dc 

11 

bdb 

28 

da 

1 

cb 

1 

dc 

2 

dd 

9 

cb 

15 

ab 

5 

ad 

9 

cd 

30 

cb 

11 

dc 

14 

ada 

14 

bb 

5 

da 

5 

dd 

6 

be 

19 

cb 

19 

cb 

22 

dc 

4 

ad 

10 

bd 

10 

ca 

22 

da 

Ave 


380,000 
170,000 
280,000 
670,000 
240,000 
130,000 
140,000 

94,000 
130,000 

50,000 
290,000 
270,000 
260,000 

70,000 
100,000 
110,000 

38,000 
480,000 
180,000 
480,000 
140,000 
140,000 
130,000 
120,000 

212,170 
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TABLE  B-2 

Transmissivity  Values  for 
Wells  Located  on  the  Bozeman  Fan 


Township 


Range 


Section 


Location  in 
Sections 


Transmissivity 
(gal/day/ft) 


IS 
2S 
2S 
2S 
2S 
2S 


5E 

26 

5E 

11 

5E 

11 

5E 

15 

5E 

27 

5E 

35 

da 
dec 
dec 
aa 
cc 
dc 


Ave 


64,000 
4,500 
36,000 
26,000 
50,000 
65,000 
40,917 
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APPENDIX  C 
As  described  in  the  report,  groundwater  return  flow  factors  were  cal- 
culated for  the  Canyon  Ferry  system.   Table  C-1   on  the  following  page 
summarized  these  calculations. 
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TABLE  C-1 
Groundwater  Return  Flow  Factors  for  the  Canyon  Ferry  Drainage 


T 
a  -   -- 

120,000  gal/day/ft 
.20 
.9  ft  /day)  30.4  days 
X  5280  ft/mi) ^ 

1  ft^ 

^7.48  gal 

-  80,213.90  ft 

/day 

at   (80213 
17  "  (2  mi  : 

-  .0218673 

at 

Return  Flow 

Factor 

Month 

17       R 

R 

X  Month 

First  Diff. 

Second  Diff 

: . 

0 

.0 

.0 

.0 

1 

.0218673 

.2224734 

.2224734 

.2224734 

.22247 

2 

.0437346 

.3146052 

.6292104 

.4067370 

.18426 

3 

.0656019 

.3849690 

1.1549069 

.5256965 

.11896 

4 

.0874692 

.4433043 

1.7732171 

.6183102 

.09261 

5 

.1093366 

.4931402 

2.4657012 

.6924841 

.07417 

6 

.1312039 

.5363153 

3.2178921 

.7521909 

.05971 

7 

.1530712 

.5740272 

4.0181906 

.8002985 

.04811 

8 

.1749385 

.6071568 

4.8572541 

.8393635 

.03907 

9 

.1968058 

.6363950 

5.7275547 

.8703006 

.03094 

10 

.2186731 

.6623034 

6.6230335 

.8954788 

.02518 

11 

.2405404 

.6853458 

7.5388036 

.9157701 

.02029 

12 

.2624077 

.7059104 

8.4709246 

.9321210 

.01635 

13 

.2842751 

.7243250 

9.4162249 

.9453003 

.01318 

14 

.3061424 

.7408674 

10.3721434 

.9559185 

.01062 

15 

.3280097 

.7557747 

11.3366200 

.9644766 

.00856 

16 

.3498770 

.7692493 

12.3079884 

.9713684 

.00689 

17 

.3717443 

.7814656 

13.2849146 

.9769262 

.00556 

18 

.3936116 

.7925734 

14.2663206 

.9814060 

.00448 

19 

.4154789 

.8027019 

15.2513364 

.9850158 

.00361 

20 

.4373462 

.8111963 

16.2392622 

.9879258 

.00291 

21 

.4592136 

.8204540 

17.2295335 

.9902713 

.00235 

22 

.4810809 

.8282588 

18.2216925 

.9921590 

.00189 

23 

.5029482 

.8354510 

19.2153731 

.9936806 

.00152 

24 

.5248155 

.8420949 

20.2102784 

.9949053 

.00122 

25 

.5466828 

.8482469 

21.2061735 

.9958951 

.00099 

26 

.5685501 

.8539564 

22.2028654 

.9966919 

.00080 

27 

.5904174 

.8592667 

23.2001997 

.9973343 

.00064 

28 

.6122847 

.8642162 

24.1980530 

.9978533 

.00052 

29 

.6341520 

.8688387 

25.1963230 

.9982700 

.00042 

30 

.6560194 

.8731643 

26.1949292 

.9986062 

.00034 

31 

.6778867 

.8772194 

27.1938023 

.9988731 

.00027 

32 

.6997540 

.8810280 

28.1928968 

.9990945 

.00022 

33 

.7216213 

.8846111 

29.1921659 

.9992691 

.00017 

34 

.7434886 

.8879876 

30.1915781 

.9994122 

.00014 

35 

.7653559 

.8911744 

31.1911038 

.9995257 

.00011 

36 

.7872232 

.8941867 

32.1907207 

.9996159 

.00009 

Note  that  the  return  flow  factor  (the  percentage  of  a  given  month's  irrigation 
groundwater  recharge  that  returns  in  a  specified  month  decreases  with  time. 
After  36  months,  99.96  percent  of  a  given  month's  recharge  has  returned  to  the 
surface  water  system.   Return  flow  in  subsequent  months  Is  assumed  to  be  zero. 
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APPENDIX  D 

Canyon  Ferry  Operations  Record 

This  appendix  contains  the  Bureau  of  Reclamation's  monthly  reservoir 
operation  records  for  the  period  from  1955  through  198A. 
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APPENDIX  E 

Table  E-1  lists  the  month-end  storage  goals  that  were  targeted  in  the 
operation  of  Canyon  Ferry  Reservoir. 

TABLE  E-1 

Month  Month-End  Storage  Goal  (Acre-ft) 

Oct  1,830,000 

Nov  1,850,000 

Dec  1,780,000 

Jan  1,680,000 

Feb  1,600,000 

Har  1,500,000 

Apr  1,500,000 

May  1,680,000 

June  1,927,000 

July  1,883,500 

Aug  1,800,000 

Sept  1,760,000 
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APPENDIX  F 

This  appendix  contains  the  Bureau  of  Reclamation's  Canyon  Ferry  power 
plane  characteristic  curves  (and  corresponding  equations) .  It  also  contains 
information  on  reservoir  surface  area-volume-elevation  relationships  and 
reservoir  storage  allocations. 
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This  appendix  contains  the  Bureau  of  Reclamation's  general  operating 
principles  for  Canyon  Ferry, 
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GENERAL  OPERATING  PRINCIPAL?  AT  CANYON  FERRY 


1.  The  cop  3  feet  between  elevation  3797  and  3800  are  used  exclusively 
for  downstream  flood  control  and  when  storage  rises  into  this  pool 
operation  of  the  reservoir  is  directed  by  the  Corps  of  Engineers. 
This  storage  is  generally  evacuated  as  fast  as  downstream  conditions 
permit. 

2.  As  soon  as  the  storage  has  peaked,  usually  in  June  or  July,  power 
releases  are  adjusted  so  that  the  pool  will  be  drawn  to  about  elevation 
3783  by  the  following  March  1.   Inflows  are  reestimated  each  month 

and  releases  adjusted  accordingly.   Releases  to  meet  this  schedule 
are  limited  to  powerplant  capacity.   Water  is  not  spilled  to  provide 
this  drawdown. 

3.  Most  of  the  stored  water  that  will  be  released  from  Hebgen  is 
released  in  October  and  November  and  rescorage  of  this  water  in 
Canyon  Ferry  may  cause  its  pool  to  rise  slightly  in  these  months. 
However,  KPC  will  try  to  limit  the  Hebgen  drawdown  during  these  months 
in  an  effort  to  maintain  Canyon  Ferry  pool  below  elevation  3794  after 
December  1.   Storage  below  elevation  379A  prior  to  winter  freeze-up 

is  desired  to  prevent  ice-jam  problems  at  the  head  end  of  the  lake. 

4.  Beginning  about  January  5,  and  monthly  thereafter  through  June  5, 
forecasts  are  made,  from  the  snow  cover  and  precipitation,  of  the 
estimated  April-July  inflow.   When  these  forecasts  become  available, 
operational  procedures  are  changed.   Based  on  forecasts  made  about 
January  5,  February  5,  and  March  5,  releases  are  set  based  on  the 
most  probable  April-July  inflow  forecast  to  allow  the  reservoir  to 
fill  to  elevation  3797  near  the  end  of  June.   There  are  two  limitations 
on  release  rates  during  these  3  months,  (a)  release  no  more  than 
powerplant  capacity,  and  (b)  release  no  less  than  required  to  produce 
critical  period  energy  at  Canyon  Ferry  and  The  Montana  Power  Company 
(MFC)  Missouri  River  powerplants.   In  a  dry  year  item  (b)  may  result 

in  the  reservoir  pool  not  reaching  elevation  3797,  the  Cop  of  the  Joint 
use  pool. 

5.  With  forecasts  made  on  or  after  April  5,  operations  are  still  based 
on  reaching  elevation  3797.   However,  if  the  forecasts  show  that 
releases  in  excess  of  powerplant  capacity  must  be  made,  the  amount 
of  water  scheduled  for  spilling  is  based  on  the  assured  runoff  (90 
percent  chance  chat  inflow  will  exceed  this)  rather  than  the  most 
probable  runoff.   Water  forecast  as  surplus  in  April  and  May  is 
scheduled  for  release  by  June  5  and  that  forecast  as  surplus  on 
June  5  usually  scheduled  for  release  by  June  15.   Spilling  of  addi- 
tional water  is  made  only  to  the  extent  that  current  inflows  and 
reservoir  contents  indicate  that  additional  spills  are  required. 
Releases  are  limited  to  15,000  c.f.s.  or  full  downstream  channel 
capacity  as  long  as  space  is  available. 
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In  years  of  unusually  heavy  snowpack  it  may  be  desirably  to  begin 
releasing  assured  spills  prior  co  April.   If  this  assured  spill  can 
be  beneficially  used  by  MPC  during  March  and  early  April  it  may  be 
released  as  long  as  it  does  not  result  in  a  loss  of  powerplant 
capability  at  Canyon  Ferry  or  jeopardize  the  refill  of  Canyon  Ferry 
Reservoir. 

Depending  on  when  the  spring  runoff  starts,  the  release  of  water  based 
on  inflow  forecasts  may  draw  the  pool  as  low  as  elevation  3770.   In  a 
series  of  dry  years  the  pool  may  be  drawn  as  low  as  elevation  3728  to 
meet  firm  power  generation  requirements,  and  satisfy  MPC  prior  water 
rights.   If  storage  is  drawn  below  elevation  3728,  the  powerplant  will 
become  inoperable.   Below  elevation  3783  total  powerplant  generating 
capability  is  limited  by  reduced  head  on  the  power  turbines. 

The  operation  of  Canyon  Ferry  is  closely  coordinated  with  MPC.   During 
low  runoff  years  it  may  be  determined  beneficial  to  the  Bureau  and 
MPC  to  generate  more  than  critical  period  energy.   This  can  be  done 
if  both  the  Bureau  and  MPC  agree  but  it  will  be  at  the  risk  of  possibly 
having  less  than  critical  period  energy,  at  a  later  date,  if  dry 
conditions  should  persist. 
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